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Abstract 

We study linear problems Sa defined on tensor products of Hilbert spaces with 
an additional (anti-) symmetry property. We construct a linear algorithm that uses 
finitely many continuous linear functionals and show an explicit formula for its worst 
case error in terms of the eigenvalues A of the operator W\ = S\S\ of the univariate 
problem. Moreover, we show that this algorithm is optimal with respect to a wide 
class of algorithms and investigate its complexity. We clarify the influence of different 
(anti-) symmetry conditions on the complexity, compared to the classical unrestricted 
problem. In particular, for symmetric problems with Ai < 1 we give characterizations 
for polynomial tractability and strong polynomial tractability in terms of A and the 
amount of the assumed symmetry. Finally, we apply our results to the approximation 
problem of solutions of the electronic Schrodinger equation. 
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1 Introduction 



In the theory of linear operators : H d — > Gd denned between Hilbert spaces it is well-known 
that we often observe the the so-called curse of dimensionality if we deal with <i-fold tensor 
product problems. That is, the complexity of approximating the operator Sd by algorithms 
using finitely many pieces of information increases exponentially fast with the dimension d. 

*This is an extended version of a same- named paper by the author which was published in the Journal of 
Approximation Theory [8]. Here all the proofs, as well as some additional assertions, are explicitly included. 
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In the last years there have been various approaches to break this exponential dependence 
on the dimension, e.g., we can relax the error definitions. Another way to overcome the curse 
is to introduce weights in order to shrink the space of problem elements Hd- In the case of 
function spaces this approach is motivated by the assumption that we have some additional 
a priori knowledge about the importance of several (groups of) variables. 

In the present paper we describe an essentially new kind of a priori knowledge. We 
assume the problem elements / G Hd to be (anti-) symmetric. This allows us to vanquish 
the curse and obtain different types of tractability. 

The problem of approximating wave functions, e.g., solutions of the electronic Schrddinger 
equation, serves as an important example from computational chemistry and physics. In 
quantum physics wave functions \I/ describe quantum states of certain (i-particle systems. 
Formally, these functions depend on d blocks of variables yj, which represent the spacial 
coordinates and certain additional intrinsic parameters, e.g., the spin, of each particle within 
the system. Due to the Pauli principle, the only wave functions \l/ which are physically 
admissible are antisymmetric in the sense that ^(y) = (— l)' 7r '^ r (7r(y)) for all y and all 
permutations ir on a subset / C {l,...,d} of particles with the same spin. Here (— l)' 7 "' 
denotes the sign of ir. The above relation means that \& only changes its sign if we replace 
particles by each other which possess the same spin. For further details on this topic we 
refer to Section 5 of this paper and the references given there. Inspired by this application 
we illustrate our results with some simple toy examples at the end of this section. 

To this end, let Hi and Gi be infinite dimensional separable Hilbert spaces of univariate 
functions /:DcR4l and consider a compact linear operator Si : Hi — > Gi with singular 
values a = (aj)j(zfq. Further, let A = (Aj)j e ^ = (cr|) jeN denote the sequence of the squares 
of the singular values of Si. Finally, assume Sd- Hd — >■ Gd to be the <i-fold tensor product 
problem. We want to approximate Sd by linear algorithms using a finite number of continuous 
linear functionals. 

By n cnt (e, d) we denote the minimal number of information operations needed to achieve 
an approximation with worst case error at most e > on the unit ball of Hd- The integer 
n ent (e,d) is called information complexity of the entire tensor product problem. Further, 
consider the subspace of all / G Hd that are fully symmetric, i.e., 

f(x) = f(ir(x)) for all x G D d and all permutations ir of {1, . . . , d}. 

The minimal number of linear functionals needed to achieve an ^-approximation for this 
subspace is denoted by n sym (e,d). Finally, define the subspace of all functions / G H d that 
are fully antisymmetric by the condition 

f(x) = (-l) kl /(7r(x)) for all x G D d and all vr 
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and denote the information complexity with respect to this subspace by n asy (e,d). 

Since is a Hilbert space, the optimal algorithm for the entire tensor product prob- 
lem is well-known. Moreover, it is known that its worst case error, and therefore also the 
information complexity, can be expressed in terms of A, i.e. in terms of the squared singular 
values of the univariate problem operator Si, see, e.g., Sections 4.2.3 and 5.2 in Novak and 
Wozniakowski [4] . It turns out that this algorithm, applied to the (anti-) symmetric problem, 
calculates redundant pieces of information. Hence, it can not be optimal in this setting. 

In preparation for our algorithms, Section 2 is devoted to (anti-) symmetric subspaces 
in a more general fashion than in this introduction. Moreover, there we study some basic 
properties. In Section 3 we conclude formulae of algorithms for linear tensor product prob- 
lems defined on these subspaces. We show their optimality in a wide class of algorithms 
and deduce an exact expression for the n-th minimal error in terms of the squared singular 
values of Si. Theorem 1 summarizes the main results. Finally, we use this error formula to 
obtain tractability results in Section 4 and apply them to wave functions in Section 5. 

Our results yield that in any case (if we deal with the absolute error criterion) 

n asy {e, d) < n sym {e, d) < n cnt {e, d) for every e > and all deN, 

where for d = 1 the terms coincide, since then we do not claim any (anti-) symmetry. 
To see that additional (anti-) symmetry conditions may reduce the information complexity 
dramatically consider the simple case of a linear operator Si with singular values o such that 
Ai = A2 = 1 and Aj = for j > 3. Then the information complexity of the entire tensor 
product problem can be shown to be 

n ent (e, d) = 2 d for all d G N and e < 1. 

Hence, the problem suffers from the curse of dimensionality and is therefore intractable. On 
the other hand, our results show that in the fully symmetric setting we have polynomial 
tractability, because 

n sym (e, d) = d+l for all d e N and e < 1. 

It can be proved that in this case the complexity of the fully antisymmetric problem decreases 
with increasing dimension d and, finally, the problem even gets trivial. In detail, we have 

n asy (e, d) = max {3 - d, 0} for all d E N and e < 1, 

which yields strong polynomial tractability. 

Next, let us consider a more challenging problem where Ai = A 2 = . . . = A m = 1 and 
Aj = for every j > m > 2. For m = 2 this obviously coincides with the example studied 
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above, but letting m increase may tell us more about the structure of (anti-) symmetric 
tensor product problems. In this situation it is easy to check that 

n ent (e, d) = m d and n asy (5, d) = \ ^ ' d ~ ™ for every d G N and all e < 1. 

10, d > m, 

Since (™) > 2 d ~ 1 for d < [^/2j, this means that for large m the complexity in the antisym- 
metric case increases exponentially fast with d up to a certain maximum. Beyond this point 
it falls back to zero. The information complexity in the symmetric setting is much harder 
to calculate for this case. However, it can be seen that we have polynomial tractability, but 
n sym (e, d) needs to grow at least linearly with d such that the symmetric problem can not be 
strongly polynomially tractable, whereas this holds in the antisymmetric setting. The entire 
problem again suffers from the curse of dimensionality. 

The reason why antisymmetric problems are that much easier than their symmetric 
counterparts is that from the antisymmetry condition it follows that f(x) = if there exist 
coordinates j and I such that Xj = x\. Another explanation for the good tractability behavior 
of antisymmetric tensor product problems might be the initial error e™ 1 *. For every choice 
of A it tends to zero as d grows, what is not necessarily the case for the corresponding entire 
and the symmetric problem, respectively. In fact, we have 

d 

= 4t ym = Af 2 , whereas e^^X}' 2 . 

j=i 

For a last illustrative example consider the case Ai = 1 and Xj + i = j" 13 for some j3 > 
and all j G N. That means that we have the two largest singular values o\ = 02 of S\ 
equal to one. The remaining series decays like the inverse of some polynomial. If /3 — the 
operator Si is not compact, since (A m ) m6 pj does not tend to zero. Hence, all the information 
complexities are infinite in this case. For /3 > 0, any 5 > and some C > it is 

n cnt (e, d) > 2 d , n sym (e, d)>d+l and n asy (e, d) < Ce- (2/fS+5) , for all e < 1, d G N. 

Thus, again for the entire problem we have the curse, whereas the antisymmetric problem 
is strongly polynomially tractable. Once more, the symmetric problem can shown to be 
polynomially tractable. Note that in this example the antisymmetric case is not trivial, 
because all Xj are strictly positive. If we replace j" 13 by log" 1 ^' + 1) in this example we 
obtain (polynomial) intractability even in the antisymmetric setting. 

Altogether these examples show that exploiting an a priori knowledge about (anti-) sym- 
metries of the given tensor product problem can help to obtain tractability, but it does not 
make the problem trivial in general. We conclude the introduction with a partial summary 
of our main complexity results. 
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Theorem. Let A = (\ m ) m eN denote the non-increasing sequence of the squared singular 
values of Si : Hi — > G\ and assume A2 > 0. Then for the information complexity of (anti-) 
symmetric linear tensor product problems Sd we obtain the following characterizations: 

• The fully symmetric problem is strongly polynomially tractable w.r.t. the normalized 
error criterion iff A G l T for some r > and Ai > A2. Furthermore, in the case Ai < 1 
the problem is strongly polynomially tractable w.r.t. the absolute error criterion iff 
A G £ T and A 2 < 1. 

• The fully antisymmetric problem is strongly polynomially tractable w.r.t. the absolute 
error criterion iff A G i T for some r > 0. 

In contrast, it is known, see Novak and Wozniakowski [4], that 

• the entire tensor product problem is never (strongly) polynomially tractable w.r.t. to 
normalized error criterion. Moreover, the problem is strongly polynomially tractable 
w.r.t. the absolute error criterion iff A G i T for some r > and Ai < 1. 

2 Spaces with (anti-) symmetry conditions 

Motivated by the example of wave functions in Section 1, we mainly deal with function 
spaces in this section. To this end, we start by defining (anti-) symmetry properties for 
functions which will lead us to orthogonal projections, mapping the function space onto its 
subspace of (anti-) symmetric functions. It will turn out that these projections applied to 
a given basis in the tensor product Hilbert function space lead us to handsome formulae 
for orthonormal bases of the subspaces. In a final remark we generalize our approach and 
define (anti-) symmetry conditions for arbitrary tensor product Hilbert spaces based on the 
deduced results for function spaces. 

We use a general approach to (anti-) symmetric functions, as it can be found in Section 2.5 
of Hamaekers [1]. Therefore, for a moment, consider an abstract separable Hilbert space F 
of real-valued functions defined on a domain Q C M. d . In this part of the paper let d > 2 be 
fixed. The inner product on F is denoted by (•, ■) F . Moreover, let / = 1(d) C {1, . . . , d} be 
an arbitrary given non-empty subset of coordinates. Then we define the set 

Si = {it : {1, . . . , d} — > {1, . . . , d} \ ir bijective and tt\^ = id} 

of all permutations on {l,...,d} that leave the complement of / fixed. Obviously, the 
cardinality of this set is given by #<S/ = (#/)!, where # denotes the number of elements of 
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a set. For a given n £ Sj we define the mapping 



7r': n ->■ M d , a; = (x ± , ...,x d )^ n'(x) = . . . ,2;^)). 



To abbreviate the notation we identify n and 7r' with each other. 

For an appropriate definition of partial (anti-) symmetry of functions / G F we need the 
following simple assumptions. For every n G Si we assume 

(Al) x G f2 implies 7r(x) G f2, 

(A2) f E F implies /(tt(-)) G F and 

(A3) there exists c n > (independent of /) such that ||/(7r(-)) | F\\ < c n \\f \ F\\. 

Note that these assumptions always hold if F is a d-fold tensor product Hilbert space H d = 
Hi <S> ■ ■ ■ <S> Hi equipped with a cross norm, as described in the examples of the previous 
section. 

Now we call a function / G F partially symmetric with respect to I (or I-symmetric for 
short) if a permutation n G Si applied to the argument x does not affect the value of /. 
Hence, 



Moreover, we call a function / G F partially antisymmetric with respect to I (or I- 
antisymmetric, respectively) if / changes its sign by exchanging the variables Xi and Xj with 
each other, where i,j G /. That is, we have 



where \ir\ denotes the inversion number of the permutation it. The term (— l)' 71 "' therefore 
coincides with the sign, or parity of n and is equal to the determinant of the associated 
permutation matrix. In the case = 1 we do not claim any (anti-) symmetry, since the 
set Si = {id} is trivial. For I — {1, ... : d} functions / which satisfy (1) or (2), respectively, 
are called fully ( anti- ) symmetric. 

Note that, in particular, formula (2) yields that the value f(x) of (partially) antisym- 
metric functions / equals zero if Xi = Xj with i 7^ j and i,j G /. For (partially) symmetric 
functions such an implication does not hold. Therefore, the (partial) antisymmetry property 
is a somewhat more restrictive condition than the (partial) symmetry property with respect 
to the same subset /. As we will see in the next sections this will also affect our complexity 
estimates. 



f(x) = f(ir(x)) for all x G f2 and every n G Si. 



(1) 



f( x ) = (-l)^f(n{x)) for all x G Q and every tt G Si, 



(2) 
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Next, we define the so-called symmetrizer & f and antisymmetrizer 2lf on F with respect 
to the subset I by 



Gf:F-^F, f^ef(f) 



1 



£/(*(•)) 



#5/ 



and 



2tf:F^F, /^2tf(/) 



1 



#5/ 



If there is no danger of confusion we use the notation (5/ and 21/ instead of ©f and 2lf , 
respectively. The following lemma collects together some basic properties which can be 
proved easily. For details see the appendix of this paper. 

Lemma 1. Both the mappings Pi G {©/,2l/} define bounded linear operators on F with 
Pj = Pi. Thus, 6/ and 21/ provide orthogonal projections of F onto the closed linear 
subspaces 



6j(F) = {f e F\f satisfies (1)} and 2t/(F) = {/ G F | / satisfies (2)} (3) 



Note that the notion of partially (anti-) symmetric functions can be extended to more 
than one subset J. Therefore, consider two non-empty subsets of coordinates I, J G {1, . . . , d} 
with I fl J = 0. Then we call a function / G F multiple partially (anti-) symmetric with 
respect to I and J if / satisfies (1), or (2), respectively, for / and J. Since I and J are 
disjoint we observe that n o a = a o n for all 7r G iS/ and a E Sj. Thus, the linear projections 
P/ G {©/, 21/} and Pj G {6j, 2lj} commute on P, i.e. P/ o Pj = Pj o P 7 . 

Further extensions to more than two disjoint subsets of coordinates are possible. We 
will restrict ourselves to the case of at most two (anti-) symmetry conditions, because in 
particular wave functions can be modeled as functions which are antisymmetric with respect 
to / and and J = I c , where I c denotes the complement of / in {1, ... , d}; see, e.g., Section 5 
of this paper. 

Up to this point the function space P was an arbitrary separable Hilbert space of (i-variate 
real-valued functions. Indeed, for the definition of (anti-) symmetry we did not claim any 
product structure. On the other hand, it is also motivated by applications to consider tensor 



of all partially (anti-) symmetric functions (w.r.t. /) in P, respectively. Hence, 



F = 6/(P) © (6/(P)) x = 2t/(P) © (2l/(P)) 



(4) 
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product function spaces; see, e.g., Section 3.6 in Yserentant [10]. In detail, it is well-known 
that so-called spaces of dominated mixed smoothness, e.g. W / 2 1 '"' ,1 ^(R 3d ), can be represented 
as certain tensor products; see Section 1.4.2 in Hansen [2]. 

Nevertheless, if we take into account such a structure, i.e., assume F = H d = Hi®. . .<g>Hi 
(d times), where Hi is a suitable Hilbert space of functions / : D — > R, it is known that we 
can construct an orthonormal basis (ONB) of F out of a given ONB of Hi. In fact, if 
{r]i | i G N} is an ONB of the underlying Hilbert function space Hi then the set of all d-fold 
tensor products {rj djj = (g)f =1 rj^ | j = . . . , j d ) G N d }, 

d 

Vd,j(x) = Y[vji(xi), x = (xi, ...,x d )e D d , 
1=1 

is mutually orthonormal in H d and forms a basis. To exploit this representation we start 
with a simple observation. 

Let j G N d and x G D d , as well as a non-empty subset / of {1, . . . , d} be arbitrarily fixed. 
If we define a = n^ 1 G Sj then 

(%imj)(x) = ^T(-i)^W^)) = -L ]T(-i) w n % m few) 

= ipr EHl^^iW = ipr £(-i) |CT ~Va)(*) (5) 
# s i ^ ™=i # s * ^ 

= -^E(- 1 ) k W(,-)W- 

For simplicity, once again we identified = n(ji, . . . ,j d ) with • • • ,j n (d)) for multi- 

indices j G N d . Obviously, the same calculation can be made for (5/ without the factor (—1). 
Since x G D d was arbitrary we obtain 

GiVdj = itV Yl ^) and a ^' = wE(- 1 ) H ^0) fora11 ^ nd - ( 6 ) 

Note that in general, i.e. for arbitrary j G N d and a G Si, the tensor products 77d,er(j) an< ^ %j 
do not coincide, because taking the tensor product is not commutative in general. Therefore, 
&i is not simply the identity on {r] d j \ j G N d }. On the other hand, we see that for different 
j G N d many of the functions &ir) d j coincide. Of course the same holds true for %lir] d j, at 
least up to a factor of (—1). 
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We will see in the following that for Pj G {&i, 21/} a linearly independent subset of all 
projections {Pirjdj \j G N d } equipped with suitable normalizing constants can be used as an 
ONB of the linear subspace Pj(Hd) of J-(anti-) symmetric functions in Hd- To this end, we 
need a further definition. For fixed d > 2 and I C {1, . . . , d}, let us introduce a function 

M I = M Itd :N d ^{0,...,#I}* 1 

which counts how often different integers occur in a given multi-index j G N d among the 
subset / of coordinates, ordered with respect to their rate. To give an example let d — 7 
and I = {1,...,6}. Then M Ij7 applied to j = (12,4,4,12,6,4,4) G N 7 gives the #7 = 6 
dimensional vector Mij(j) = (3, 2, 1, 0, 0, 0), because j contains the number "4" three times 
among the coordinates ji, . . . "12" two times and so on. Since in this example there are 
only three different numbers involved, the fourth to sixth coordinates of Mjj(j) equal zero. 
Obviously, Mj is invariant under all permutations 7r G Si of the argument. Thus, 

M 7 (j) = M/(7r(j)) for all j G N d and tt G 5/. 

In addition, since Mi(j) is again a multi-index, we see that \Mi(j)\ = j^l and Mi{j)\ are 
well-defined for every j G N d . With this tool we are ready to state the following assertion 
which can be shown using elementary arguments as well as Lemma 1; see the appendix. 

Lemma 2. Assume {r/dj \j G N d } to be a given orthonormal tensor product basis of the 
function space H d and let ^ I — . . . , i#j} C {1, . . . , d}. Moreover, for Pj G {©/, 21/} 
define functions : JJ a! — > R, 



Then the set | fc G V^} builds an orthonormal basis of the partially (anti-) symmetric 
subspace Pi(Hd), where Vd is given by 

{k eN d \k h < k i2 < . . . < k t#1 }, if Pj = Sj, 
{k G N d \ k h < k l2 < ... < k i#I }, ifP 7 = Stj. 

Observe that in the antisymmetric case the definition of ^ for j G simplifies, since 
then M/(j)! = 1 for all j G V^. Moreover, note that in the special case I — {1, . . . , #7} we 
have 
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That is, we can consider the subspace of /-(anti-)symmetric functions / G H d as the tensor 
product of the set of all fully (anti-) symmetric ^/-variate functions with the (d — #/)-fold 
tensor product of H\. Modifications in connection with multiple partially (anti-) symmetric 
functions are obvious. 

Finally, note that Lemma 2 also holds if the index set N of the univariate basis {rji | % G N} 
is replaced by a more general countable set equipped with a total order. But let us shortly 
focus on another generalization of the previous results. 

Remark 1 (Arbitrary tensor product Hilbert spaces). Up to now we exclusively dealt with 
Hilbert function spaces. However, the proofs of Lemma 1 and Lemma 2 yield that there are 
only a few key arguments in connection with (anti-) symmetry such that we do not need this 
restriction. 

Starting from the very beginning we need to adapt the definition of I- (anti-) symmetry 
due to (1) and (2). Of course it is sufficient to define this property at first only for basis 
elements. Therefore, if E d = {rj d:k I k G ~N d } denotes a tensor product ONB of H d and 
^ I C {1, . . . , d} is given then we call an element r\ d ^ = ®f=i Vk partially symmetric with 
respect to I (or I -symmetric for short), if 

Vd,k = Vdrtk) for a11 n e S ii 

where Si and n(k) = (k n m, . . . , k n u)) are defined as above. Analogously we define I- 
antisymmetry with an additional factor (—1)'^'. Moreover, an arbitrary element in is 
called /-(anti-)symmetric if in its basis expansion every element with non- vanishing coeffi- 
cient possesses this property. 

Next, the antisymmetrizer %i is defined as the uniquely defined continuous extension of 
the linear mapping 

21,: E d -)• H d , MVd,k) = J^-i^W) 

from E d to H d . Again the symmetrizer &i is given in a similar way. Hence, in the general 
setting we define the mappings using formula (6), which we derived for the special case. Note 
that the triangle inequality yields ||P/|| < 1, for Pj G {(5/,2l/}. 

Once more we denote the sets of all J- (anti-) symmetric elements of H d by Pi(H d ), where 
Pj G {©/,2t/}. Observe, that this can be justified since the operators Pi again provide 
orthogonal projections onto closed linear subspaces as described in Lemma 1. 

Finally, also the proof of Lemma 2 can be adapted to the general Hilbert space case. 
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3 Optimal algorithms 



In the present section we conclude optimal algorithms for linear problems defined on (anti-) 
symmetric subsets of tensor product Hilbert spaces as described in the previous paragraph. 
Moreover, we deduce formulae for the n-th minimal errors of these (anti-) symmetric prob- 
lems and recover the known assertions for the entire tensor product problem. 

3.1 Basic definitions and the main result 

Throughout the whole section we use the following notation. Let Hi be a (infinite dimen- 
sional) separable Hilbert space with inner product (-,-) H and let G± be some arbitrary 
Hilbert space. Furthermore, assume Si : Hi — > Gi to be a compact linear operator between 
these spaces and consider its singular value decomposition. That is, define the compact 
self-adjoint operator W\ = S\Si : H ± — > Hi and denote its eigenpairs with respect to a 
non-increasing ordering of the eigenvalues by {(e^, Aj) | i G N}, i.e. 

Wi(ei) = Xiei, and (ei,ej) Hi = 5 id with Ai > A 2 > . . . > 0. (8) 

Then A = (Aj) i6 N coincides with the sequence of the squared singular values a 2 = (of ) ieN 
of Si and the set {e; | % G N} forms an ONB of Hi, see, e.g., Section 4.2.3 in Novak and 
Wozniakowski [4]. In the following we will refer to Si as the univariate problem or univariate 
case. 

For d > 2, let H d — Hi <g> . . . <g> Hi be the ci-fold tensor product space of H x . This means 
that H d is the closure of the set of all linear combinations of formal objects / = (S)f=i fi with 
fi e Hi, called simple tensors or pure tensors. Here the closure is taken with respect to the 
inner product in H d which is defined such that 

Id d \ d 

=tl(fi,9i) Hl for fudi^Hi. 

\ 1=1 1=1 I H d 1=1 

With these definitions H d is also an infinite dimensional Hilbert space and it is easy to check 
that 

Ed = = (g) Vn e H d | j = (ji, . . . , d) e N d | (9) 

forms an orthonormal basis in H d if {r]i e Hi \ % e N} is an arbitrary ONB in the underlying 
space Hi. Similarly, let G d = Gi <E> • • • <S> Gi, d times, and define S d as the tensor product 
operator 

S d = Si ® . . . ® Si : H d G d . 
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In detail, we define the bounded linear operator S d : E d — > G d such that for all j G N d we 
have Sd(i] d j) = ^((S^zLi ^jj = 0f =1 5'i(^j ; ) G Gd- Then Sd is assumed to be the uniquely 
defined linear, continuous extension of Sd from Ed to H d . 

We refer to the problem of approximating Sd- Hd — > G d as the entire d-variate problem. 
In contrast, we are interested in the restriction S d \ , Pi(H d ) — > Gd of S d to some (anti-) 
symmetric subspace Pi(H d ) with P I G as described in the previous section. To 

abbreviate the notation we denote this restriction again by S d and refer to it as the I- 
( anti-) symmetric problem. 

For the singular value decomposition of the entire problem operator S d we consider the 
self-adjoint, compact operator 

W d = S d ] S d : H d ^H d . 

Its eigenpairs {(e d j, Xdj) \j — C?i, • • ■ ,jd) are given by the set of all ci-fold (tensor) 

products of the univariate eigenpairs (8) of W 1: i.e., 

d d 

£d,j = (g) e n and \ dJ = JJ \ n for j = . . . , j d ) G N d . (10) 
1=1 i=i 

It is well-known how these eigenpairs can be used to construct a linear algorithm A' n d which 
is optimal for the entire <i-variate tensor product problem. In detail, A! d minimizes the 

worst case error 

e™(A n4 ;H d )= sup \\A n , d (f) ~ S d (f) \ G d \\ 
feB(H d ) 

among all adaptive linear algorithms A n>d using n continuous linear functionals. Here B{H d ) 
denotes the unit ball of H d . In other words, A' n d achieves the n-th minimal error 

e(n,d;H d )= inf e WOT (A M ; H d ). 

With this notation our main result reads as follows. 

Theorem 1. Let {(e m , X m ) m G N} denote the eigenpairs of W\ given by (8). Moreover, for 
d > 1 let ^ / = {ii, . . . , i#r} C {1, . . . , d} and assume S d to be the linear tensor product 
problem restricted to the /-(anti-)symmetric subspace Pj{H d ), where Pj G {<5/,2l/}, of the 
d-io\d tensor product space H d . Finally, let be given by (7) and define 

{(^p(v),XdMv)) \veN} = {(^, A d>fc ) | k G V d } (11) 
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by £ k = y /#S I /M I (k)\-P I (e kl ®. -.®e kd ) and \ d<k = Y\ti X k n for k G V d , where ij>: N ->• V d 
provides a non- increasing rearrangement of {A^fc | G V^}. 

Then for every d > 1 the set (11) denotes the eigenpairs of W d \ p ,„ = S^Sd- Thus, for 
every n G No, the linear algorithm A* nd : Pi(H d ) — > P I (Gd), 

n 

A *n,df = ^2(f'^(v)) Hd • S dC^(v), (12) 

which uses n linear functionals, is ra-th optimal for S d on Pj{H d ) with respect to the worst 
case setting. Furthermore, it is 

e(n, d; P,(^)) = e wor (A; d ; P,(^)) = yfx^y (13) 

Let us add some remarks on this theorem. First of all, the sum over an empty index set is 
to be interpreted as zero such that A* Q d f = 0. Further, note that the worst case error can be 
attained with the function ^( n +i)- It can be improved neither by non-linear algorithms using 
continuous information, nor by linear algorithms using adaptive information. Moreover, 
observe that the classical entire tensor product problem is included as the case = 1, 
where we do not claim any (anti-) symmetry. Then = N d and the ^'s simply equal the 
tensor products e d)k = <8>f =1 e kr Hence, A* nd = A' n d . 

The remainder of this section is devoted to the proof of the main result Theorem I. 

3.2 Proof of Theorem 1 

We start with an auxiliary result which shows that any optimal algorithm A* for S d needs 
to preserve the (anti-) symmetry properties of its domain of definition, i.e. A* f G Pi(G d ) for 
all / G Pj(H d ). The following proposition generalizes Lemma 10.2 in Zeiser [11] where this 
assertion was shown for the approximation problem, i.e. for S d = id. A comprehensive proof 
can be found in the appendix of this paper. 

Proposition 1. Let d > 1 and assume ^ I C {1, . . . , d}. Furthermore, for X G {H, G}, 
let Pf denote the (anti-) symmetrizer Pi G {&i, 21/} on X d with respect to I and suppose 
A: Pf(H d ) — > G d to be an arbitrary algorithm for S d . Then, for g G H d: 

(S d oP?)(g) = (PfoS d )(g), (14) 

and for all / G P?(H d ) it holds 

\\S d f - Af | G d \\ 2 = \\S d f - Pf{Af) | G d \\ 2 + \\Af - Pf{Af) \ G d \\ 2 . (15) 
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Hence, an optimal algorithm A* for preserves (anti-) symmetry, i.e. 

A*fePf{G d ) for all f e Pf(H d ). 

Beside this qualitative assertion we are interested in explicit error bounds. Therefore, 
the next proposition shows an upper bound on the worst case error of the algorithm A* ld 
given by (12). 

Proposition 2 (Upper bound). Under the assumptions of Theorem 1 the worst case error 
of A* n d given by (12) is bounded from above by 

e™{Al d -Pj{H d ))< Jx^y 
Proof. By Lemma 2 we have for all / G Pi(H d ) the unique representation 

kev d 

Therefore, the boundedness of S d together with (11) implies that 

s d f = J2 </> &>h„ ■ S& = E </» ^w>h„ • s&w for ever y / e P '( H *)- ( 16 ) 

Furthermore, in the case Pj = 21/ it is easy to see that we have 

{S d ^S d ^ k ) Gd = - f==^ == ■ (SMje^S^e,,,)^ 

= ; 1 = V (-1)^+^1 (S d e dn(i) ,S d e da(k) ) n 

for i,j G Vd, because of the commutativity of S d and 21/ due to (14) in Proposition 1. 
Obviously, the same calculation can be done in the symmetric case, where Pj = &j. Since 
ed,ir(j) an d e d ,a(k) are orthonormal eigenelements of W d = S d ^S d : H d — > H d , see (10), it 
is \S d e d ^(j), Sd e d,a(k)) G = *d,n(j)( e dMj)i e d,<T(k)) Hd = \i,nV)8«(j),*(k)- Hence, similar to the 
proof of the mutual orthonormality of | k G V^} for Lemma 2 we obtain 

(S d £j, S d £ k ) Gd = X d ,jS jtk for all j, k G V d . (17) 

Therefore, we calculate for n G N and / G Pi(H d ) 

2 



v>n 
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On the other hand, for / G B(P I (H d )), we have by Parseval's identity 



1>||/|P,(^)I| 2 = II/I^I 



Thus, because of the non- increasing ordering of {Xd,ip(v))vem due to the choice of the rear- 
rangement tp, we can estimate the worst case error 



e™(A* n4 ; P^H,)) 2 = sup \\S d f - A* n J \ G d \\ < A^ (n+1) , 

feB(P z (H d )) 



as claimed. 



Note that formula (17) in the proof of Proposition 2 together with Lemma 2 yields that 
the set (11) describes the eigenpairs of the self-adjoint operator 

W *\ Pd H d ) = S ^ S d- P iW -+ P AHd) 

as stated in Theorem 1. Therefore, the upper bound given in Proposition 2 is sharp and 
A* hd in (12) is n-th optimal, due to the general theory; see, e.g., Corollary 4.12 in Novak 
and Wozniakowski [4]. From the general theory it also follows that adaption does not help 
to improve this n-th minimal error, see [4, Theorem 4.5], and that linear algorithms are best 
possible; see [4, Theorem 4.8]. Hence, the proof of Theorem 1 is complete. 

Since it seems to be a little bit unsatisfying to refer to these deep results for the proof 
of such an easy theorem we refer the reader to the appendix where a nearly self-contained 
proof of the remaining facts can be found. Moreover, there we describe what we mean by 
adaption in this context. 



4 Complexity 

In this part of the paper we investigate tractability properties of approximating the linear 
tensor product operator S d on certain (anti-) symmetric subsets Pj(H d ) = P Id (H d ), where 
P G {&, 21} and ^ I d C {1, . . . , d}. Therefore, as usual, we express the n-th minimal error 
derived in formula (13) in terms of the information complexity, i.e. the minimal number of 
information operations needed to achieve an error smaller than a given e > 0, 

n(s, d; Pi{H d )) = min {n G N | e(n, d; Pi{H d )) < e} . 

To abbreviate the notation we write n cnt (e,d) if we deal with the entire tensor product 
problem. Furthermore, as in the introduction, we denote the information complexity of the 
fully (anti-) symmetric problem by n asy (e,oQ and n sym (e,d), respectively. 
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4.1 Preliminaries 



From Theorem 1 we obtain for any e > and every d G N 

n(e, d; = min {n G N | A diV , (n+1) < e 2 } = # | fc G V d | JJ A fei > e 2 

by solving (13) for ip. Using this expression we can easily conclude the results for the first 
two problems in the introduction. There we dealt with the case Ai = . . . = \ m = 1 and 
Xj = for j > m > 2. 

Let us recall some common notions of tractability. If for a given problem the information 
complexity n(e, d) increases exponentially in the dimension d we say the problem suffers 
from the curse of dimensionality. That is, there exist constants c > and C > 1 such that 
for at least one e > we have 

n(e, d)>c-C d 

for infinitely many d G N. More generally, if the information complexity depends exponen- 
tially on d or e' 1 we call the problem intractable. Since there are many ways to measure the 
lack of exponential dependence we distinguish between different types of tractability. The 
most important type is polynomial tractability. We say that the problem is polynomially 
tractable if there exist constants C,p > 0, as well as q > 0, such that 

n(e, d)<C-e- p - d q for all d G N, e G (0, 1]. 

If this inequality holds with q — 0, the problem is called strongly polynomially tractable. If 
polynomial tractability does not hold we say the problem is polynomially intractable. For 
more specific definitions and relations between these and other classes of tractability see, 
e.g., the monographs of Novak and Wozniakowski [4, 5, 6]. 

In the following we distinguish two cases. First we consider the absolute error criterion, 
where we investigate the dependence of n(e,d; Pj(H d )) on \ je and on the dimension d for 
every e G (0, 1] and d G N. Note that without loss of generality we can restrict ourselves to 
e < min {l, e™ 1 *} since obviously n(e, d; Pi(H d )) = for all e > Here 

£ - = e( o, * Pl(Hd)) - ;; ; : °- 

describes the initial error of the <i-variate problem on the subspace Pi (Ha) where ip: N — > 
Vd again is a non-increasing rearrangement of the set of eigenvalues {A^ | k G V^} of 
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w d 



Pi(H d ) 



Sd'Sa and b d = d — a d denotes the number of coordinates without (anti-) 



symmetry conditions in dimension d, i.e. a d = #ld and b d = d — #I d , respectively. 

Afterwards, we deal with the normalized error criterion, where we especially investigate 
the dependence of n(e' ■ £ l d ll ,d; Pi(H d )) on \/e' for e' G (0, 1). That is, we search for the 
minimal number of information operations needed to improve the initial error by a factor e' 
less than one. 

To avoid triviality we will assume e™ 1 * > 0, for every d G N, in both cases, because 
otherwise we have strong polynomial tractability by default. From this assumption it follows 
that A x > 0, which simply means that S d is not the zero operator. Moreover, note that in 
the case of antisymmetric problems, if the number of antisymmetric coordinates, i.e. the set 
/ = 1(d), grows with the dimension, the condition e™ 1 * > (for every d G N) even implies 



Finally, we always assume A2 > 0, because otherwise S d is equivalent to a continuous lin- 
ear functional which can be solved exactly with one information operation; see Novak and 
Wozniakowski [4, p. 176]. 

For the study of tractability for the absolute error criterion we use a slightly modified 
version of Theorem 5.1, [4]. It deals with the more general situation of arbitrary compact 
linear operators between Hilbert spaces. In contrast to Novak and Wozniakowski we drop the 
(hidden) condition e™ 1 * = 1 for the initial error in dimension d. For the sake of completeness a 
proof can be found in the appendix. If we denote Riemann's zeta function by ( the assertion 
reads as follows. 

Proposition 3. Consider a family of compact linear operators {T d : F d — > G d \ d G N} be- 
tween Hilbert spaces and the absolute error criterion in the worst case setting. Furthermore, 
for d G N let (\d,i)im denote the non-negative sequence of eigenvalues of T d 'T d w.r.t. a 
non-increasing ordering. 

• If {T d } is polynomially tractable with the constants C, p > and q > then for all 
t > p/2 we have 



that 



Ai > A 2 > . . . > 0. 




l/r 



< OO, 



(18) 



where r = 2q/p and / : N ->■ N with f(d) = \(l + C) d q ] . 
In this case C T < C 2 ' p C(2r/p) 1 / T . 
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If (18) is satisfied for some parameters r > 0, r > and a function /: N — > N such 



that f(d) = 
is polynomial 



C ■ (min {e™ 1 *, l}) P ■ d q , where C > and p,q > 0, then the problem 
ly tractable and n(e, d) < (C + C T T ) r max < p - 2r } d max { 9 - rT > for every d E N 



and any e G (0, 1]. 

Let us add some comments on this result. Since, clearly, 

1 < (min {4 nit , 1}) ~ P for all p > 

Proposition 3 provides a characterization for (strong) polynomial tractability, similar to [4, 
Theorem 5.1]. But, compared with the assertions from the authors of [4], our result yields 
the essential advantage that the given estimates incorporate the initial error e™ 1 *. Hence, if 
ef 11 is sufficiently small then we can conclude polynomial tractability while ignoring a larger 
set of eigenvalues in the summation (18). 

Observe that the first statement does not cover any assertion about the initial error, since 
f(d) > 2. Hence, it might happen that we have (strong) polynomial tractability though the 
largest eigenvalue A^i = {£ l d lt ) 2 tends faster to infinity than any polynomial. To this end, 
for d € N, consider the sequences (Xd, m )meN given by 

\ d l = e 2d and \ dm = — , for m > 2. 

m 

Here, obviously, the initial error grows exponentially fast to infinity, but nevertheless the 
second point of Proposition 3 shows that {5^} is strongly polynomially tractable, since (18) 
holds with r = p = q = 0, and C = r = 2. 

Let us now return to our /-(anti-)symmetric tensor product problems S d as defined in 
Section 3. Therefore, let ^ l d — {ii, . . . , i#i d } C {1, . . . , d} and Pj d e {&i d , 2l/ d } for every 
d > 1. We start by using Proposition 3 to conclude a simple necessary condition for (strong) 
polynomial tractability of {S d } in the worst case setting w.r.t. the absolute error criterion. 
Recall that ip: N — > V d defines a rearrangement of the parameter set given in (7). That 
is, 

{\w(v) \ veN}= |A d , fe = J[ \ kl | k e V d J (19) 
denotes the set of eigenvalues of S d ( S d with respect to a non-increasing ordering, see Theorem 1. 
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Lemma 3 (General necessary conditions). The fact that {S d } is polynomially tractable with 
the constants C, p > and q > implies that A = (X m ) m em for all r > p/2. Moreover, 
for any such r and all d G N the following estimate holds: 

Proof. From Proposition 3 we know that for r > p/2 and r = 2g/p it is 

SU P^ ( E A ^) ) <0 °' ( 2 °) 



den 



,v=f(d) 



where the function / : N — > N is given by f(d) = |~(1 + C) d q ] . 

Note that for the proof of the first assertion we only need to consider the case where all 
A m are strictly positive. Then the condition (20), in particular, implies that the sum in the 
brackets converges for every fixed d G N. If we denote the subset of indices j G of the 
f(d) — 1 largest eigenvalues \ d ^( v ) by L d then there exists a natural number s — s(d) > d 
such that L d is completely contained in the cube 

Q d>a = {l,...,s} d . (21) 

Hence, we can crudely estimate the sum from below by 

oo 

E X h - E X h = E A <W) < 00 • 
j£Vd\Q d ,s jeV d \L d v=f(d) 

Since R dtS = {j = (1, 2, . . . , d — 1, to) G N d | to > s} is a subset of \ independently 
of the concrete (anti-) symmetrizer P Id , where P G {©,21}, we obtain 

oo 

(Ai • A 2 • . . . • \d-iT E y m = E A ^' ^ E A ^r 

m =s+i ieR d , s ieV d \Q djS 

Thus, for each fixed d G N the tail series ^m=s(0+i * s finite, which is only possible if 
|| A | £ T || < oo. Hence, A G £ T is necessary for (strong) polynomial tractability. 

Let us turn to the second assertion. Obviously, (20) implies the existence of some constant 
C\ > such that 

oo 

E A W) < C ^ T for a11 rf e N - 

v=f(d) 
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Indeed, Proposition 3 yields that we can take Ci = C 2r/p C(2r/p). The rest of the sum can 
also be bounded easily for any rfeN, 

f(d)-l 

due to the ordering provided by ip. Since ^ fceVd A^ = S^li ^dv(u)' ^ remams to show that 
f(d) — 1 < (1 + C*)^ for every d G N with A^(i) > 0, which is also obvious due to the 
definition of /. ■ 

Since we know that antisymmetric problems are easier than symmetric problems we have 
to distinguish these cases in order to conclude sharp conditions for tractability. 



4.2 Tractability of symmetric problems (absolute error) 

Beside the general assertion A G £ T , we start with necessary conditions for (strong) poly- 
nomial tractability in the symmetric setting. By bd we denote the amount of coordinates 
without symmetry conditions in dimension d. 

Proposition 4 (Necessary conditions, symmetric case). Let {Sd} be the problem considered 
in Lemma 3 and assume P = &. 

• If {S d } is polynomially tractable and Ai > 1 then b d G 0(\nd). 

• If {Sd} is strongly polynomially tractable and A x > 1 then b d G 0{1) and A 2 < 1/Ai. 

Proof. Assume Ai > 1 and let r be given by Lemma 3. Then, independent of the amount of 
symmetry conditions, we have \d,ip(i) — Af > 1 and there exist absolute constants r > and 
C > 1 such that 

^E A i^ C ^ deN > ( 22 ) 
Al kev d 

due to Lemma 3. In the case of strong polynomial tractability we even have r = 0. For 
d > 2 we use the product structure of Xd,k, k G Vd, and split the sum on the left w.r.t. the 
coordinates with and without symmetry conditions. Hence, we conclude 

E K**> d = a d + b d >2, 

m 1 / h£N a d, 
h 1 <...<ha d h 1 <...<h ad 

(23) 
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which leads to 



gPTd, 1=1 \ 1 / 
Z...<ha d 

In any case the second sum in the above inequality is bounded from below by 1. Thus, 
we conclude that (1 + A5/A[) 6d < (X^m=i Kn/^i) bd needs to be polynomially bounded from 
above. Since we always assume A 2 > this leads to the claimed bounds on bd in the case of 
(strong) polynomial tractability. 

It remains to show that A 2 < 1/Ai is necessary for strong polynomial tractability. To 
this end, assume for a moment A 2 > 1/Ai. Then it is easy to see that (independent of the 
number of symmetry conditions) there are at least 1 + |_cZ / 2 J different k G such that 
Ad,fe > 1. Namely, for / = 0, . . . , [d/2\ we can take the first d — I coordinates of k G Vd equal 
to one. To the remaining coordinates we assign the value two. 

In other words, we have XdMi+\d/2\) > 1- On the other hand, strong polynomial tractabil- 
ity implies X^lri+ci KiMv) — ^i f° r some absolute constants t,C,Ci > and all d G N; 
see (20). Hence, for every d > 2 [1 + C] , 

oo l+|d/2j 

Ci> E A ^w> E A ^)>A^ (1+Lci /2j)(2+^/2j-ri+ci)>^/2j + i-rci, 

v=\l+C] v=\l+C] 

because of the ordering provided by ip. Obviously, this is a contradiction. Thus, we have 
A 2 < 1/Ai and the proof is complete. ■ 

Note in passing that the previous argument can also be used to show that (independent 
of the number of symmetry conditions) the information complexity n(e, d) needs to grow at 
least linearly in d if we assume A 2 > 1/Ai. In particular, we cannot have strong polynomial 
tractability if Ai = A 2 = 1. 

We continue the analysis of /-symmetric problems with respect to the absolute error 
criterion by proving that the stated necessary conditions are also sufficient for (strong) 
polynomial tractability. To this end, we need a rather technical preliminary lemma that can 
be proven by elementary induction arguments. For the convenience of the reader we included 
also this proof in the appendix. 
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Lemma 4. Let (/i m ) me N be a non-increasing sequence of non-negative real numbers with 
fii > 0. Then, for all VgN and every d G N, it holds 



J2 »d,k<»id V 



( \ 

d 



keN d , 
i<ki<...<k d 



L=l 



(24) 



V V+2< j [V<...< j p ) 



Now the sufficient conditions read as follows. Once again, we denote the number of 
coordinates without symmetry conditions in dimension d by bd- 

Proposition 5 (Sufficient conditions, symmetric case). Let {Sd} be the problem considered 
in Lemma 3, assume P = & and let A = (A m ) m£ N G ^r f° r some r G (0, oo). 

• If Ai < 1 then {S d } is strongly polynomially tractable. 

• If Ai = 1 > A 2 and b d G 0(1) then {Sd} is strongly polynomially tractable. 

• If Ai = 1 and bd G O(lnd) then {S^} is polynomially tractable. 

Proof. Step 1. We start the proof by exploiting the property A G i T0 . It is easy to see that 
the ordering of ( y \ m )men implies 

oo 

m\Z < X? + . . . + XZ < A? = IIA I U\ T ° < oo 

i=i 

for any m G N. Hence, there exists some C To > such that A m is bounded from above by 
C To ■ vrT r for every r < 1/to- Therefore, there is some index such that for every larger m G N 
we have A m < 1. We denote the smallest of these indices by m . Similar to the calculation 
of Novak and Wozniakowski in [4, p. 180] this leads to 

OO ^OO £j T i 

J2 XT m<(P+ ^0 + C; / X ^ rdx = (P + ^0 + ^ 
Jmn+v ' ' 



m=ni 7m 0+P 



rr — 1 (mo + p) rr 1 

for every p G No and all r such that rr > 1. Thus, in particular, with r = 1/t we conclude 

1/t / r l/(l/ro-l/r)\ r(l/T0-l/r) 

E^^(p+^ + i7^(^+ir) fora " T >™ eN °- 

m=mo \ / 
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Note that for a given 5 > there exists some constant T\ > r such that for all r > r± it 
is l/(l/r — 1/t) G (r ,r + 5). Hence, if p G N is sufficiently large then we obtain for all 

r > Ti 

in=mo x 

£(p + 1)A; , (i + ^(^)*° +5> , 

ri \m +pj 

where we set Ci = max {l, Cy° +<5 }. Finally, since A mo < 1, both the summands tend to zero 
as r approaches infinity. In particular, there need to exist some r > T\ > r such that 



oo 1 

£ K. < 5. 



m=rrto 



S'tep £ All the stated assertions can be seen using the second point of Proposition 3. 
Indeed, for polynomial tractability, it is sufficient to show that 

X lk < Cd rT for all deN (25) 

feev d 

and some C, r > as well as some r > 0. If this even holds for r = we obtain strong 
polynomial tractability. 

In the case Ai < 1 we can estimate the sum on the left of (25) from above by (X)m=i Kn) d - 
Using Step 1 with m = 1 we conclude X]feev d Ki,k — 2~ d for some large r > r . Hence, the 
problem is strongly polynomially tractable in this case. 

For the proof of the remaining points assume Ai = 1. In any case XlfcgVi ^ifc — 
X^m=i Kn = II^KrolT < 00 f° r & H r ^ r 0) because of A G £ T0 . Therefore, we can as- 
sume d > 2 in the following. Again we split the sum in (25) with respect to the coordinates 
with and without symmetry conditions, i.e., for d = ad + bd > 2 we have 

(OO \ b d 

1 + A ™ I S X ^d,h- ( 26 ) 
m=2 / h£N a d, 
hi<...<ha d hi<...<h ad 

If A2 < 1 and bd is universally bounded then the first factor can be bounded by a constant 
and the second factor can be estimated using Lemma 4 with V — 0, d replaced by ad and ji 
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replaced by A T . It follows that if r is large enough we have 

a d a d / oo \ L oo 

E ^.»<i+E E ^™<i+E(E^ <i+Y,2- l = 2, 

h£N a d, L=l j( L )m L L=1 \m=2 / L=l 

hl<...<ha d 2< j [ L) <...<j[ L) 

where we again used Step 1 and the properties of geometric series. Thus, J2kev d Ki,k ls 
universally bounded in this case and therefore the problem is strongly polynomially tractable. 

To prove the last point we argue in the same manner. Now bd G 0(\nd) yields that the 
first factor in the splitting (26) is polynomially bounded in d. For the second factor we again 
apply Lemma 4, but in this case we set V = m — 2, where mo denotes the first index m G N 
such that X m < 1. Keep in mind that this index is at least two because of Ai = 1. On the 
other hand it needs to be finite, since A G £ T0 . Therefore, the second factor in the splitting 
(26) is also polynomially bounded in d due to the same arguments as above. All in all, this 
proves (25) and the problem is polynomially tractable in this case. ■ 

We summarize the results obtained for /-symmetric tensor product problems in the fol- 
lowing theorem. 

Theorem 2 (Tractability of symmetric problems, absolute error). Assume S±: Hi — > G\ 
to be a compact linear operator between two Hilbert spaces and let A = {X m ) m eN denote 
the sequence of non-negative eigenvalues of W\ = S{Si w.r.t. a non- increasing ordering. 
Moreover, for d > 1 let ^ Id C {1, . . . , d}. Assume Sj to be the linear tensor product 
problem restricted to the /^-symmetric subspace &i d (Hd) of the rf-fold tensor product space 
Hd, consider the worst case setting w.r.t. the absolute error criterion and let Ai < 1. 
Then the problem is strongly polynomially tractable if and only if A G £ T for some r > and 

• Ai < 1, or 

• 1 = A x > A 2 and (d - #I d ) G 0(1). 

Moreover, the problem is polynomially tractable if and only if A G £ T for some r > and 

• Ai < 1, or 

• Ai = 1 and (d - #I d ) G 0(\nd). 



24 



4.3 Tractability of symmetric problems (normalized error) 

Here we briefly focus on the normalized error criterion for the /-symmetric setting. Since 
(e d mt ) 2 = \d,ip(i) = for any kind of symmetric problem, this means that we have to 
investigate the influence of d and 1/e' on 

n(e' ■ ef\ d; &i d {H d )) = min {n G N \ \ d ^ n +i) < (e') 2 ^(i)} 




for e' G (0,1), de N. 



Hence, in fact we have to study the information complexity of a scaled tensor product problem 
S' d : &i d {H d ) — > G d with respect to the absolute error criterion. The squared singular values 
of S[ equal \i = (// m ) m gN with fj, m = X m /Xi. Obviously, we always have /ii = 1. Furthermore, 
/i G £ T if and only if A G l T . This leads to the following theorem. 

Theorem 3 (Tractability of symmetric problems, normalized error). Assume S\ \ Hi — >■ G 1 
to be a compact linear operator between two Hilbert spaces and let A = (A m ) me N denote 
the sequence of non-negative eigenvalues of W\ = S\Si w.r.t. a non- increasing ordering. 
Moreover, for d > 1 let ^ Id C {1, . . - ,d}. Assume S d to be the linear tensor product 
problem restricted to the /^-symmetric subspace &i d (H d ) of the <i-fold tensor product space 
H d and consider the worst case setting w.r.t. the normalized error criterion. 
Then the problem is strongly polynomially tractable if and only if A G i T for some r > 
and Ai > A 2 and (d - #/ d ) G 0(1). 

Moreover, {Sd} is polynomially tractable if and only if A G £ T for some r > and {d — ^I d ) G 
0(h\d). 

4.4 Tractability of antisymmetric problems (absolute error) 

We start this subsection with simple sufficient conditions for strong polynomial tractability. 

Proposition 6 (Sufficient conditions, antisymmetric case). Let {Sd} be the problem con- 
sidered in Lemma 3, assume P = 21 and let A = (A m ) me N G £ TQ for some r G (0, oo). 

• If Ai < 1 then {S d } is strongly polynomially tractable, independent of the number of 
antisymmetry conditions. 

• If Ai > 1 and if there exist constants r > r and d G N such that for the number of 
antisymmetric coordinates a d in dimension d it holds that 

hi^l > lnfllAKX) for all d > d (27) 
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then the problem {S d } is also strongly polynomially tractable. 

Proof. Like for the symmetric setting, the proof of these sufficient conditions is based on the 
second point of Proposition 3. We show that under the given assumptions 

E X dMv) < C < cx) for every d G N 

yV = l J 

and some r > t . Once again ip and Vd are given as in (19) and (7), respectively. 
Since for d = 1 there is no antisymmetry condition we have ip = id and 





E A i-^) = E A * = iiakii<iiak i 



Therefore, due to the hypothesis A G £ To the term for d = 1 is finite. 

Hence, let d > 2 be arbitrarily fixed. For s G N with s > d we define the cubes of 
multi-indices similar to (21). With this notation we obtain the representation 



oo 



Yl ^Mv) - Yl x ik ~ J™ Yl x ik- 

v=i kev d fceV d nQ diS 

Without loss of generality we may reorder the set of coordinates such that Id — . . . ,i ad } — 
{1, . . . , a<f}. That is, we assume partial antisymmetry with respect to the first ad coordinates. 
Furthermore, we define U ad;S = {j G Q ad , s \ji<j 2 <---< ja d } and set b d = d - a d . 

If bd > then the set of multi-indices under consideration splits into two non-trivial 
parts: 

V d D Q d ^ s = U ad;S x Q bdtS for all s > d. 
Because of the product structure of A^fc (k G V ' d) this implies 

(ad \ ( b d 

e e m 
3&Ua d ,s 1 = 1 ) \i£Qb d ,s 1 = 1 

Since the sequence A = (A m ) meN is an element of t T we can easily estimate the second factor 
for every s > d from above by 

b d b d s / s \ b d / oo \ yr-b d -r 

E U X l=UY X m= E A - <[Y XT m) =HA|4ir. (28) 

1=1 m=l \m=l / \m=l / 
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To handle the first term we need an additional argument. Note that the structure of U adjS 
implies 

a d a d a d 

E II E IK+<^ ! E II 

j&Qa d ,s 1=1 3&Qa d ,s 1 = 1 j^ U a d ,s 1=1 

3k,m:j k =j m 

which leads to the upper bound 

O-d -i a d 

E IK<^ E n^iiAKir-, 

j£Ua d , s 1 = 1 a ' j€Qa d ,s 1 = 1 

where we used the same arguments as in (28). Once again this upper bound does not depend 
on s > d. Hence, due to d = ad + b d , we conclude 

E A w,=,'i- E ^<^nAK T ir 

v=i feev d nQ d , s 

for every choice of 2l/ d . Of course, for every 2 < d < do this upper bound is trivially less 
than an absolute constant. Thus, we can assume d > d . Then, due to the hypothesis of the 
second point we have ln^!) > In (||A | ^11^), which implies 

E A W)J < (^H A I £ r\\ Td ) <1 for d>d . 

Hence, (27) is sufficient for strong polynomial tractability, independent of Ai. Therefore it 
suffices to show that Ai < 1 implies (27) in order to complete the proof. To this end, let 
Ai < 1. We know from Step 1 in the proof of Proposition 5 that there exists some r > t 
such that 



OO 

iiai4-ii t = XX<-<i. 



m=l 



Thus, we see that the right hand side of (27) is negative, whereas the left hand side is 
non- negative for every choice of a^. ■ 

We also briefly comment on this result. First, note that a sequence A = (A m ) mg N that 
is not included in any £ r -space, < r < oo, has to converge to zero more slowly than the 
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inverse of any polynomial, i.e., m~ a for a > arbitrarily fixed. Thus, only sequences like 
A m = l/ln(m) lead to polynomial intractability in the fully antisymmetric setting. 
Secondly, observe that (27) is quite a weak assumption. For example if we have 



a d > 



d 



In d a 



with < a < 



ln(||AK r |D 



for all sufficiently large d then 

Mgj) > aJMaJ, - I) > 1 In (X ^) _^ l > ^ ^ ^ 

a a a ma a 

If ct equals its upper bound, i.e. a = l/ln(||A | ^-ID, then the condition (27) does not hold. 
This also shows that assumptions like ad = \d^ with j3 < 1 are not sufficient to conclude 
(27). 

Note that Proposition 3 allows us to omit the largest f(d) — 1 eigenvalues \d,tp( v ) where 
f(d) may grow polynomially in (e™*) -1 with d. We did not use this fact in the proof of the 
sufficient conditions. 

Let us now turn to the necessary conditions. We will see that we need a condition 
similar to (27) in order to conclude polynomial tractability if we deal with slowly decreasing 
eigenvalues A. 

Proposition 7 (Necessary conditions, antisymmetric case). Let {Sd} be the problem con- 
sidered in Lemma 3 and assume P = 21. Furthermore, let {Sd} be polynomially tractable 
with the constants C, p > and q > 0. 

Then, for d tending to infinity, the initial error ef^ tends to zero faster than the inverse of 
any polynomial. Furthermore, A = (A m ) mg N £ for every r > p/2 and for all 5 > there 
exists some do G N such that 



i / ha i f ir\ 

m(||All r |D -8 < iVln ' for all d > d . 

d k=i \ X k J 



(29) 



Thus, we have Ai < 1 or lim^oo ad = oo. 

Proof. Step 1. For the whole proof assume r > p/2 to be fixed. Then Lemma 3 shows that 
A G i T . Moreover, we again use the notation d = ad + bd, where ad = #Id denotes the number 
of coordinates with antisymmetry conditions in dimension d. Similar to the symmetric case 
we can split the sum of the eigenvalues such that for all d e N 

/ OO \ b d 

E ^ fc = E ^n) E ^d,j - 11^ I ^r\\ Tl>d • A[ • . . . • A„ 



fcev d \m=i / jeN a d, 

l<jl<-<ja d 
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Hence, Lemma 3 implies that 



W^ + >* + ^(7)(S)' 

In what follows we will use this inequality to conclude all the stated assertions. 

Step 2. Here we prove the limit property for the initial error e™ 1 * = a/A^i), i.e. we need 
to show that for every fixed polynomial V > 

\ d ^i)V(d) ^ if d ->■ oo. (31) 

Since A^i) = A^ d • Ai • . . . X ad < X\ d ■ X\ d = Xf due to the ordering of A = (X m )meN we can 
restrict ourselves to the non-trivial case Ai > 1 in the following. Now assume that there 
exists a subsequence (dkjken of natural numbers as well as some constant Cq > such that 
Xd k ,^(i)V(dk) is bounded from below by Co for every fceN. Then for every d = dk the right 
hand side of (30) is bounded from above by some polynomial V\{dk) > 0. On the other 
hand, due to the general condition A 2 > 0, the term ||A | £ T \\ T /X\ is strictly larger than one. 
Thus, it follows that there exists some C\ > such that 

bd k < C\ ln(dfc) for every k e N. 

Therefore, we have in particular ad k = dk — bd k — > oo, for k — > oo. Moreover, the assumed 
boundedness of Xd k ^(i)V(d k ) leads to 



CoVidk)' 1 < X dkMl) < X? ln{dk) - Ax •.. . X a , = d c k lHXl) • Ai • . . . A 



since Ai > 1. As we showed in Step 1 of the proof of Proposition 5 the fact A G i T yields 
the existence of some C T > such that X m < C T m~ l l T for every m G N. Indeed, this 
holds for C T = ||A|£ T ||, which needs to be larger than one because of Ai > I. Hence, 
A[ • . . . • Xl dk < C T T adk (a dfc !) _1 , what implies 

(^f) adk <a dk \<(c;r^v 2 (d k ), ken 

for some other polynomial V2 > 0. Thus, if k is sufficiently large we conclude 

a dk <a dk ]n(j^j < ln(P 2 (4)), 

since a dk — > 00 implies a dk /(eC^) > e for k > ko. Therefore, the number of antisymmetric 
coordinates a d needs to be logarithmically bounded from above for every d out of the sequence 
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(dk)k>k - Because also bd k was found to be logarithmically bounded this is a contradiction 
to the fact d\~ = a& k + h^ k . Thus, the hypothesis \d k ^(\)P{dk) > Co > can not be true for 
any subsequence (dk)k- In other words it holds (31). 

Step 3. Next we show (29). From the former step we know that there needs to exist some 
d* E N such that l/\d,tp(i) > 1 f° r & h d> d*. Hence, (30) together with r > p/2 implies 

M*r\r\*^(**v = c 2d ^ for ^ 



where we set C 2 = 1 + C + C 2r / p ((2r/p). Therefore, we conclude 

\ \\- d < ft l|A|4l|T 

2 fc=i fc 

for all d> d* , which is equivalent to 

hiww-^4D"(^). (32) 



k=l 



Obviously, for given 5 > 0, there is some d** such that ln(C 2 d 2qT/p )/d < 5 for all d > d**. 
Hence, we can choose d = max {d*, d**} in order to obtain (29). 

Step 4- It remains to show that Ai > 1 implies that lim^oo ad is infinite. To this end, 
note that every summand in (32) is strictly positive. If we assume for a moment the existence 
of a subsequence (rffc)feeN such that ad k is bounded for every fceN then the right hand side of 
(32) is less than some positive constant divided by dk- Hence, it tends to zero if k approaches 
infinity. On the other hand, for large d, the left hand side of (32) is strictly larger than some 
positive constant, because of Ai > 1 and A 2 > 0. This contradiction completes the proof. ■ 

As mentioned before there are examples such that the sufficient condition (27) from 
Proposition 6 is also necessary (up to some constant factor) in order to conclude polynomial 
tractability in the antisymmetric setting. Now we are ready to give such an example. 

Example 1. Consider the situation of Lemma 3 for P = 21 and assume the problem {Sd} 
to be polynomially tractable. In addition, for a fixed r G (0, oo), let A = (A m ) mS N £ be 
given such that Ai > 1 and, moreover, assume that there exist m G N such that 

A m > - — for all m > m and some a > 1. (33) 



m 
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Then we claim that for every 5 > there exists rfeN such that 



(^-S^ln{\\X\e T \\ T ) < for all d > d. (34) 

Recall that due to Proposition 6, for the amount of antisymmetry ad, it was sufficient to 
assume 

In (|| A | £ T \\ T ) < ^ n ^\^' for every d larger than some fixed d G N 
a 

in order to conclude strong polynomial tractability. 

Before we prove the claim it might be useful to give a concrete example where (33) holds 
true. To this end, set \ m = 1/m 2 , r = 1, a = 3 and m = 2. Then it is easy to check that 
11-^ I ^t\\ = C(2) = tt 2 /6 and obviously we have Ai = 1. 

To see that the claimed inequality (34) holds true we can use Proposition 7 and, in 
particular, inequality (32). Since Ai > 1 we know that lim^a^ = oo, i.e. > m for every d 
larger than some d\ G N. Moreover, note that (33) is equivalent to 

In (- — | T ^ ^ < adn(m) for all m > mo- 



Hence, if d > d 1 we can estimate the right hand side of (32) from above by 

if, f\\M^r\\ T \^rno /||A|4|r\ a f , , C x hx(a d \) 

- d E >» (-xj-J £ t ■ ln ) * , > , lnW £ t + a — ■ 

k=l \ K s \ m / k=m Q +l 

Consequently, this leads to 

^niM^lin Ca + ln(C 2 d^/P) ^ ln(q d !) 

— ln( A £ r ) ; < ; — for d > max \d , d\\ . 

a a ■ d d 

Now (34) follows easily by choosing d > maxjfi*,^} large enough such that the negative 
term on the left is smaller than a given 5 > times ln(||A | £ T 



Although there remains a small gap between the necessary and the sufficient conditions 
for the absolute error criterion, the most important cases of antisymmetric tensor product 
problems are covered by our results. We summarize the main facts in the next theorem. 

Theorem 4 (Tractability of antisymmetric problems, absolute error). Let Si : Hi — > G\ be a 
compact linear operator between two Hilbert spaces and let A = (A m ) me ^ denote the sequence 
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of non-negative eigenvalues of W\ = S\Si w.r.t. a non-increasing ordering. Moreover, for 
d > 1 let 7^ 1^ C {1, . . . , d}. Assume S d to be the linear tensor product problem restricted 
to the /^-antisymmetric subspace 2t/ d (H d ) of the ci-fold tensor product space H d and consider 
the worst case setting with respect to the absolute error criterion. 
Then for the case Ai < 1 the following statements are equivalent: 

• {S d } is strongly polynomially tractable. 

• {Sd} is polynomially tractable. 

• There exists a universal constant r e (0, oo) such that A G £ T . 

Moreover, the same equivalences hold true if Ai > 1 and j^l d grows linearly with the dimen- 
sion d. 

Finally, before we continue with the normalized error criterion, we want to deduce an 
exact formula for the complexity in the case of fully antisymmetric functions. Hence, we set 
I — I d — { 1 , . . . , d} for every dGN and consider 

in the following. In this case the set of parameters is given by 

V d = {k = (h, . . . , k d ) e N d | h < k 2 < . . . < k d }. 

Thus, to obtain a worst case error less or equal than a given e > we need at least 

n asy (5, d) = n(e, d; %(H d )) = # j k G V d | X d , k = J[ \ h > 5 2 | 

linear functionals for the rf-variate case. 

Due to the ordering of (X m ) me ^ the largest eigenvalue A^^ with k e is given by the 
square of e™ 1 * = ■ X 2 ■ . . . ■ X d - In other words, it is 

n as >-(e,rf) = for all e > e™*. 

To calculate the cardinality also for e < e 1 ^ let us define 

i d {5 2 ) = min [i e N | a { = X { ■ X i+1 ■ . . . ■ X i+d -i < 5 2 } , for 5 > and deN. (35) 

Using this notation we can formulate the following assertion. Again the proof can be found 
in the appendix of this paper. 
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Proposition 8 (Complexity of fully antisymmetric problems). Let {S^} be the problem 

considered in Lemma 3 and assume P = 21 as well as / = — {1, . . . , d}. 

Then for every e > the information complexity is given by n asy (e, d) = n ont (e, 1), if d — 1, 

and 

n asy (e,rf) (36) 

i d (e 2 ) id-i(£ 2 A ;i -i) i2(e 2 /[\ h - 1 -..,Xi d _ 2 - 1 ]) 

= E E - E 

il=2 ?2='l + l id_l=id_2+l 

if d > 2. Here the quantities ij, for j = 2, . . . , d, are defined as in (35). 
Remark 2. If we define 

fc-i 

Z=0 

for fc G N and a non-increasing sequence Ai > A2 > . . . > 0, then we can interpret the 
quantities ik{5 2 ) as information complexities of modified univariate problems S[ k \ In detail, 
let S[® : H\ — > Gi define a compact linear operator such that 

Wi k) = (S[ k) y (S[ k) ) : H ± ^ H! 

possesses the eigenvalues {at' | m G N}. Then 

n cnt (5,l)=n cnt (5,l;S { 1 h) : H 1 ^G 1 )=i k (5 2 )-l for all 5 > 0. 

Further, note that for k > 2 the quantities i^{e 2 /{Xi^-i ■ . . . • A/ d A ._i]) are non-increasing 
functions in l ± , . . . , and e. 

Out of Proposition 8 we can conclude bounds on the information complexity. If d > 2 
and £ < e™ 1 * then the sum in (36) contains at least the term with the index li — 2, l 2 — 
3, . . . , ld-i = d. That is, for any choice of A we get the lower bound 

n asy (e, d) > n cnt (e/y/X 1 - ... ■ A d _i, l) - d + 1, 

which can be used to show that we cannot expect the same nice conditions for (strong) 
polynomial tractability as before if we switch from the absolute to the normalized error 
criterion. We conclude this subsection with a corresponding example. 



n cnt (e/ V / A Zl _ 1 -...-A id _ 1 _ 1 , l) - l d ^ + 1 
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Example 2. Assume X m = m 2a for all m G N and some a > 0. Then we need to estimate 

n ent (e/y/\ 1 -...-\ d . 1 , l) = n cnt (e ■ (d - l)\ a , 1) 

= # {m G N | m- 2a >(e-(d- l)\ a ) 2 } 

= #<mGN|m< —. \- -1 > —. — - 1. 

Therefore, 

™ asy ( £ > d ) > -T, h — ttt - d = d ( - - 1 ] if d > 2 and £ < e™*. 

v y ~ eVa ■ (d- 1)! V £ / °- rf! / ~ 

Since in this case the initial error e™ 1 * for the ci-variate problem equals l/d\ a , we need at 
least 

n"V • 4 nit ,rf) > d - l) 

linear functionals to improve the initial error by a factor e' < 1. Because this bound grows 
linearly with the dimension the problem is not strongly polynomially tractable with respect to 
the normalized error criterion. Nevertheless, the sequence A is an element of l\/ a , say, which 
implies strong polynomial tractability for the absolute error criterion due to Theorem 4. 

4.5 Tractability of antisymmetric problems (normalized error) 

Up to now every complexity assertion in this paper was mainly based on Proposition 3 
which dealt with the general situation of arbitrary compact linear operators between Hilbert 
spaces and with the absolute error criterion. While investigating tractability properties of 
/-symmetric problems with respect to the normalized error criterion, we were able to use 
assertions from the absolute error setting. Since for /-antisymmetric problems the struc- 
ture of the initial error is more complicated, this approach will not work again. Therefore, 
we start this subsection with a modified version of another known theorem by Novak and 
Wozniakowski [4, Theorem 5.2]. 

Proposition 9. Consider a family of compact linear operators {T^: — > | d G N} 
between Hilbert spaces and the normalized error criterion in the worst case setting. Further- 
more, for d G N let (\d,i)ien denote the non-negative sequence of eigenvalues of T d ^T d w.r.t. 
a non-increasing ordering. 
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• If {T d } is polynomially tractable with the constants C,p > and q > then for all 
r > p/2 we have 



where r = 2q/p and /: N — > N with f(d) = 1. 
In this case, C£ < 1 + C + C 2r/p C(2r/p). 

• If (37) is satisfied for some parameters r > 0, r > and a function / : N — > N such that 
/(d) = \Cd q ~\, where C > and g > 0, then the problem is polynomially tractable 
and n(e' ■ e^, d) < (C + C T T ) (e')~ 2r d m ^ q ' rT \ 

Note that this shows that (strong) polynomial tractability is characterized by the bound- 
edness of the sum over the normalized eigenvalues, were we are allowed to omit the Cd q 
largest of them. Of course, our results are equivalent to the assertions given by Novak and 
Wozniakowski [4], as one can see easily. But now the connection between the different error 
criterions is more obvious. From this point of view Proposition 9 reads more natural than 
[4, Theorem 5.2]. The key is to apply the same proof technique for both the assertions. 

Moreover, observe that also the theorem in [4] for the normalized error criterion includes 
further assertions concerning, e.g., the exponent of strong polynomial tractability. Again 
our proof implies the same results. 

Similar to the former sections we continue with an application of Proposition 9 to our 
antisymmetric tensor product problems. To this end, assume Si : H 1 — > G\ to be a compact 
linear operator between two Hilbert spaces and let A = (\ m )m&i denote the sequence of 
non-negative eigenvalues of W\ = S\Si w.r.t. a non-increasing ordering. Moreover, for d > 1 
let 7^ Id — {1, • • • , d}. Assume Sd to be the linear tensor product problem restricted to the 
/^-antisymmetric subspace %i d {H d ) of the rf-fold tensor product space H d and consider the 
worst case setting w.r.t. the normalized error criterion. Finally, let hd denote the number 
of coordinates without antisymmetry conditions in dimension d, i.e. bd = d — ad, where 
a d = #I d for d e N. 

Proposition 10 (Necessary conditions, antisymmetric case). Under these assumptions the 
fact that {Sd} is polynomially tractable with the constants C, p > and q > implies that 
A E i T for all r > p/2. 

Moreover, for d tending to infinity, £™ rt tends to zero faster than the inverse of any polynomial 
and bd € 0(\nd). Thus, lim^oo a,d/d = 1. 

In addition, if {Sd} is strongly polynomially tractable then b d G 0(1). 




(37) 
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Proof. From Proposition 9 it follows that there is some C\ > such that for every d G N 



d,m fcev rf v=i V A ^(i)/ 



if r > p/2. Once more the rearrangement function ^ and the index set are given as in (19). 
Indeed, the proof of Proposition 9 yields that it is sufficient to take C\ = l+C+C 2T / p ((2r/p). 
In particular, for d — 1 it is Vi = N and Ai^ = A&, for k G N, such that we have ip = id 
because of the ordering of A. Hence, we conclude 

oo 

\\X\£ T \\ T = < CiX[ < oo. 

In other words, A G l T . Moreover, like with the arguments of Step 1 in the proof of 
Proposition 7, it follows 



(38) 



since \d,ip(i) — A? d • Ai • . . . • \ ad and ||A | l T \\ T > \{ due to the general assertion A 2 > 0. Thus, 
polynomial tractability of {Sj} implies bd < C 2 ln(rf) for some C 2 > 0, i.e. bd G 0(lnd). 
Therefore, obviously, we have 

1 -^ a d , b d In d Ind 

1 > — - = 1 - — • — — > 1 - C 2 ■ — > 1, d ->■ oo. 

a Ind d d 

The proof that strong polynomial tractability leads to b d G (9(1) can be obtained using 
(38) with the same arguments as before and q = 0. Finally, we need to show the assertion 
concerning e™ 1 *. To this end, we refer to Step 2 in the proof of Proposition 7. ■ 



5 Application: wave functions 

During the last decades there has been considerable interest in finding approximations of 
wave functions, e.g., solutions of the electronic Schrodinger equation. Due to the so-called 
Pauli principle of quantum physics only functions with certain (anti-) symmetry properties 
are of physical interest. 

In this last section of the present paper we briefly introduce wave functions and show 
how our results allow to handle the approximation problem for such classes of functions. For 
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a more detailed view, see, e.g, Hamaekers [1], Yserentant [10], or Zeiser [11]. Furthermore, 
for a comprehensive introduction to the topic, as well as a historical survey, we refer the 
reader to Hunziker and Sigal [3] and Reed and Simon [7]. 

In particular, the notion of multiple partial antisymmetry with respect to two sets of 
coordinates is useful for describing wave functions In computational chemistry such func- 
tions occur as models which describe quantum states of certain physical rf-particle systems. 
Formally, these functions depend on d blocks of variables yi = (x^\ s^), for i = 1, . . . , d, 
which represent the spacial coordinates \ ) £ M. and certain additional 

intrinsic parameters £ C of each particle y within the system. Hence, rearranging the 
arguments such that x = (x^\ . . . , x^) and s = (s^\ . . . , s^) yields that 

^: (R 3 ) d x C d ->■ R, (x,s) V(x,s). 

In the case of systems of electrons one of the most important parameters is called spin and 
it can take only two values, i.e., E C = { — Due to the Pauli principle the only 
wavefunctions \1/ that are physically admissible are those which are antisymmetric in the 
sense that for / C {1, . . . , d} and I c = {1, . . . , d} \ I 

*(7r(x),7r(s))) = (-l) W ^(:r,s) for all neS I US I c. 

Thus, ^ changes its sign if we replace any particles yi and yj by each other which posses 
the same spin, i.e. = So, the set of particles, and therefore also the set of spacial 
coordinates, naturally split into two groups i+ and /„. In detail, for wave functions of d 
particles yi we can (without loss of generality) assume that the first #/ + indices % belong 
to the group of positive spin, whereas the rest of them possess negative spin, i.e. I + = 
{l,...,#/ + }and/_ = {#/ + + l,...,4. 

In physics it is well-known that some problems, e.g., the electronic Schrodinger equation, 
which involve (general) wave functions can be reduced to a bunch of similar problems, where 
each of them only acts on functions fy s out of a certain Hilbert space F d = F d (s). That is, 

V s = V(., s )eF d = {f: (R 3 ) d ^R} 

with a given fixed spin configuration s G C d . Of course, every possible spin configuration s 
corresponds to exactly one choice I + C {1, . . .,d} of indices. Moreover, it is known that 
F d is a Hilbert space which possesses a tensor product structure. Therefore, we can model 
wave functions as elements of certain classes of smoothness, e.g., F d C H d = W / ^ 1 '"'' 1 ' ) (1R M ), 
as Yserentant [10] recently did, and incorporate spin properties by using the projections of 
the type 21 = o 2l/_, as defined in Section 2. 
In particular, Lemma 2 yields 

F d = K{H d ) = %i + (H #I+ ) <g) % Z _(H #I _) 
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and the system of all 

lk = y/#Si + -#S I _-*(r lk ), keW d , 

with 

V d = {A; = G \ h<i 2 < ... < i#i + and ji < . . . < j #I _} 

builds an orthonormal basis of Fd = Ql(Hd), where the set {7% | k — (ki, . . . , k d ) G N d } is once 
again assumed to be an orthonormal tensor product basis of Hd — Hi <g) . . . <g) H\ constructed 
with the help of {rji \ i G N}, an arbitrary orthonormal basis of Hi. 

Note that in the former sections the underlying Hilbert space H 1 always consists of 
univariate functions. In contrast wave functions of one particle depend on three variables, 
but we want to stress the point that this is just a formal issue. However, this approach 
radically decreases the degrees of freedom and improves the solvability of certain problems 
Sd like the approximation problem, i.e. Si = id: Hi — > G±, considered in connection with 
the electronic Schrodinger equation. 

Theorem 1 then provides an algorithm which is optimal for the G^- approximation of 
d-particle wave functions in Fd with respect to all linear algorithms that use at most n 
continuous linear functionals. Moreover, the error can be calculated exactly in terms of the 
squared singular values A = (A m ) m£ M of Si. 

Furthermore, it is possible to prove a modification of Theorem 4 for problems dealing with 
wave functions. In fact, for the mentioned approximation problem polynomial tractability 
as well as strong polynomial tractability are equivalent to the fact that the sequence A of 
the squared singular values of the univariate problem belong to some £ T -space if we consider 
the absolute error criterion. The reason is that all the assertions in Section 4.4 can be easily 
extended to the multiple partially antisymmetric case. In detail, if we denote the number 
of antisymmetric coordinates within each antisymmetry group I™ C {1, . . . , d} by a d:m , 
m = 1, . . . , M, then the constraint ad + bd = d extends to a d ,i + . . . + a d ,M + b d = d. Here bd 
again denotes the number of coordinates without any antisymmetry condition. In conclusion, 
the sufficient condition (27) in Proposition 6 transfers to 

1 M 

-^2\n(a d;m \) >\\X\£ T \\ T , for all d > d , 

m=l 

which is always satisfied in the case of wave functions, since then M = 2 and at least 
the cardinality a d m of one of the groups of the same spin needs to grow linearly with the 
dimension d. 
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Appendix 



Proof of Lemma 1 in Section 2 

We show Lemma 1 in the case of function spaces. 

Proof. Obviously, P G {6/,2l/} is well-defined due to the assumptions (Al) and (A2). The 
linearity directly follows from the definition and, using (A3), the operator norm is bounded 
by max {c n \ n G Si}. 

To show that the operators are idempotent, i.e. P 2 = P, we first prove that 2l/(/) 
satisfies (2) for every / G F. Therefore, we use the representation 

(*/(/))(*(*)) = £(-l) H /(<7(7r(*))) = E(-l) |A|+W /(A(x)) = (-l)W(2t,(/))(x) 

0-65/ Ae«Sj 

for every fixed 7r G Sj. Here we imposed \ = a o n e Si and used 

|A o 7r _1 | = |A| + Itt" 1 ! = |A| + |tt| . 

Hence, we have 21/ (F) C {/ G F \ f satisfies (2)}. In a second step, it is easy to check that 
for every function g G F which satisfies (2) it is 21/ (p) = g. Thus, {/ G F \ f satisfies (2)} C 
2l/(F) and 21/ is a projector onto 21/ (F). 

Since the same arguments also apply for the symmetrizer (3/ this shows (3), as well as 
P 2 = P for P G {6/, 21/}. Because of the boundedness of the operators the subsets ©/(F) 
and 21/ (F) are closed linear subspaces of F and we obtain the orthogonal decompositions 

F = &i(F) © (6/(F)) x = 2t/(F) © (2l/(F)) ± , 

where the X denotes the orthogonal complement with respect to (•, -) F , i.e. the image of the 
projectors (id — ©/) and (id — 21/), respectively. ■ 

The proof of Lemma 1 in the case of arbitrary tensor product Hilbert spaces works exactly 
in the same way. 

Proof of Lemma 2 in Section 2 

We prove Lemma 2 in the case of function spaces. For the case of arbitrary tensor prod- 
uct Hilbert spaces only slight modifications are needed. Indeed, the only difference is the 
conclusion of formula (39) in Step 2. In the general setting this simply follows from our 
definitions. 
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Proof. Step 1. We start by proving orthonormality. Therefore, let us recall (6). To abbreviate 
the notation further, we suppress the index H d at the inner products (•, -) H in this proof. 
For Pi = 21/ and j, k G easy calculus yields 



Assume now there exists me{l,...,d} such that j m ^ k m . Then the ordering of j, k G V<f 
implies that ^ a(k) for all a, it G <Sj, since 7r and a leave the coordinates m E I c fixed. 
Hence, we conclude that we have ir(j) = a(k) only if j = k. 

At this point we have to distinguish the antisymmetric and the symmetric case. For 
P = 21/ the only way to conclude 7r(j) = o~(k) is to claim j = k and n = a. Furthermore, we 
see that in the antisymmetric case we have Mi(j)\ = 1 for all j E V<f, since all coordinates ji, 
where / G /, differ. Therefore, in this case the last inner product coincides with Sj^ • 5 VjtJ 
because of the mutual orthonormality of {rj d j | j G N d }. Hence, 



So, let us consider the case Pj = &i and j = k G V d , because we already saw that 
otherwise (Cj,Ck) equals zero. Then for fixed a G Sj there are M/(j)! different permutations 
7r G Sj such that = cr(j). This leads to 



and completes the proof of orthonormality. 

Step 2. It remains to show that the span of {£ fc | k G V^} is dense in Pi(H d ) for P G 
{@/,2l/}. To this end, note that every multi-index j G N d can be represented by a uniquely 
defined multi-index k G and exactly Mi{k)\ different permutations n G Si such that 
j = n(k). 

Now assume / G %li(H d ), i.e. / G H d satisfies (2). Then the expansion of f(ir(x)) with 
respect to the basis functions {r/ d j | j G N d } in H d yields for n E Si 





for all j, k G V d 



as claimed. 




(_1)M (y ) . = ^ (/, 77 d)fc ) • 77 d;fc (7r(x)) for every x G 
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Similar to the arguments used in (5) we have i] d ^(n(x)) = r] d;a ^(x) with a = n 1 . Therefore, 
we conclude for x G D d 

Yl (( _1 ) H ■ •feW = E (f,Vdrt*(k))) -Vd,a(k)(x) = (f>Vd,*V))-ridj(x). 

jGN d feGN d jGN d 

Because the expansion is uniquely defined we get 

(-1) H -{f,Vdj) = (f,VdMJ)) f or all 3 eN d and vr G S T . (39) 

Using the observations from the beginning of this step we can decompose the basis expansion 
of / e %i(H d ) C H d and use the derived formula (39) to get 

jeN d kev d aeS! n >' feev d A >' a^Si 

Now (6) yields that 

f = E y iff)! • (/ ' % ' fc) • y iff)! ' ^ (% ' fe) - 

fe^V^ V V 

Furthermore, summing up (39) with respect to 7r leads to 

(/,%,*> = ^ ^(-i) w = //,-^ £(-i) |w W>) = </, , 

7rG<Sj \ 7rG<S/ / 

for k G Vd, such that finally / G 2l/(if rf ) possesses the following representation 

/= E (/,&>•&, 

since = a/ #Si/Mi(k)\ ■ %Li(r)d,k) per definition. This proves the assertion for the case 
Pi — 21/ . The remaining case Pj = (5/ can be treated in the same way. ■ 

Proof of Proposition 1 in Section 3.2 

Proof. The proof is organized as follows. First we show that the problem operator S d and 
the (anti-) symmetrizer Pj commute on H d , i.e. it holds (14). In a second step we conclude 
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(15) out of this. The (anti-) symmetry of A*f for an optimal algorithm A* then follows 
immediately. 

Step 1. Assume E d = {f]d,j \ j G N d } to be an arbitrary tensor product ONB of Hd, as 
defined in (9). Then, for fixed j G N d , formula (6) and the structure of the linear tensor 
product operator Sd = Si ® . . . ® Si yields in the case Pi = 21/ 



\^° 7 neS! 1=1 



(0 

i=i 

= inr E(- 1 ) H <g) s ^ Ml) ) = %?(s d (v d ,)). 

^ 1 7re«Sj 1=1 

Obviously, the same is true for Pj = 6/. Hence, it holds (14) at least on the set of basis 
elements E d of H d . Because of the representation 

9=J2{g,Vd,j) Hd -r]d,j, geH d , 

as well as the linearity and boundedness of the operators Pf , Pf and Sd we can extend the 
relation (14) from Ed to the whole space Hd- 

Step 2. Now let / G P^(Hd) and let Af denote an arbitrary approximation of Sdf ■ Then 
Sdf = Sd^Pf 1 f) = Pp(Sdf), due to Step 1. Using the fact that Pf provides an orthogonal 
projection onto Pf(Gd), see (4), we obtain (15), 

\\Sdf-Af\G d \\ 2 = \\P?(S d f)-[Pf(Af) + (id G -Pf)(Af)]\G d \\ 2 

= \\P?(S d f - Af) | G d \\ 2 + ||(id G - P?){Af) | G d \\ 2 

= \\S d f - Pf(Af) | G d \\ 2 + \\Af- Pf{Af) | G d \\ 2 , 
as claimed. ■ 



A self-contained proof of Theorem 1 in Section 3 

In order to deduce a lower bound on the n-th minimal error of approximating S d on (anti-) 
symmetric subspaces Pj(H d ), where P G {21, ©}, let us define the classes of algorithms under 
consideration. 

An algorithm A n ^ for Sd'- Fd = Pi(Hd) — > Gd which uses n pieces of information is 
modeled as a mapping <fi : W 1 — > Gd and a function N : Fd — > M. n such that A n ^ = ° N. In 
detail, the information map iV is given by 

N(f) = (Li(f), L 2 (f), . . . , L n (f)) , / G F d , (40) 
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where Lj G A. Here we distinguish certain classes of information operations A. In one case 
we assume that we can compute continuous linear functionals. Then A = A al1 coincides 
with F*[, the dual space of F^. If L v depends continuously on / but is not necessarily linear 
the class is denoted by A cont . Note that in both the cases also A" is continuous and we 
obviously have A al1 C A cont . 

Furthermore, we distinguish between adaptive and non-adaptive algorithms. The latter 
case is described above in formula (40), where L v does not depend on the previously computed 
values Li(f), . . . , L„_i(/). In contrast, we also discuss algorithms of the form A n ^ = (p o N 
with 

N(f) = (Lx(/), L 2 (f; Vl ), . . . , L n (f; y u . . . , y n _0) , f E F d , 

where y\ = Li(f) and y v = L v (f] yi, . . . , y v -i) for v — 2, 3, . . . , n. UN is adaptive we restrict 
ourselves to the case where L v depends linearly on /, i.e. L v ( • ; y±, . . . , y v -i) G A al1 . 

In all cases of information maps, the mapping <fi can be chosen arbitrarily and is not 
necessarily linear or continuous. The smallest class of algorithms under consideration is the 
class of linear, non-adaptive algorithms of the form 

n 

A n ,df = • g v , 

v=l 

with some g v G and L v G A al1 . We denote the class of all such algorithms by A 1 ™- On the 
other hand, the most general classes consist of algorithms A n ^ = <fi o N, where is arbitrary 
and either uses non-adaptive continuous or adaptive linear information. We denote the 
respective classes by A™ nt and ^4 adapt . 

For the proof that the upper bound given in Proposition 2 is sharp we use a generalization 
of Lemma 1 in W. [9] . 

Lemma 5. Suppose S to be a homogeneous operator between linear normed spaces X and 
Y, i.e. S(ax) = aS(x) for all x G X and a G M.. Furthermore, assume that V C X is a 
linear subspace with dimension m and there exists a constant a > such that 

a.||/|X||<||S(/)|y|| for all / G V. 

Then for every n < m and every algorithm A n G A^ nt U ^4 adapt 

e WOI (A n ; S,X) = sup \\S(f) - A n (f) \ Y\\ > a. 
feB(x) 
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Proof. It is well-known that for A n = <fi o iV with n < m there exists /* G V such that 
N(f*) = N(-f*) and ||/ | X\\ = 1. Thus, A n (f*) = A n (-f*). For a more detailed view, see, 
W. [9, Lemma 1] and the references in there. Using the triangle inequality for Y we obtain 

e wor (A n ;X) > m a x{\\S(±n-A n (±n\Y\\} = m^{\\S(r)±A n (f*)\Y\\} 

> Imn + Mn i *i + \\s(n - Mn i *i) 

> i||2.5(r)|F||>a||r|X||=a 

and the proof is complete. ■ 

Now let X = Pji^Hd) and Y = Gd- Furthermore, for a given n G N , define a = \d,%i>{n+i) 
and consider V = span {£^(1), . . . , £^( n +i)} C Pi(Hd). Then, obviously, dimU = n + 1 = 
m > n . With the representation (16) and formula (17) from the proof of Proposition 2 we 
conclude 

n+1 n+1 

\\Sdf \G d \\ 2 = ^2(f,^( v )) • -Vv-o) > ^ (/, £v>0)) = fl2 11/ I -^ll 2 ' / e ^ 

D=l D=l 

where we used the monotonicity of {Ad^^eN and Parseval's identity. This leads to the 
desired lower bound result: 

Proposition 11 (Lower bound). Under the assumptions of Theorem 1 the n-th minimal 
error with respect to the class A^ nt U ^ dapt is bounded from below by 

e(n, d; Pj(H d )) = inf e wor (A M ; Pj{H d )) > J X dMn+ i) for all d G N, n G N . 

Hence, this together with Proposition 2 shows that A* nd given in (12) is n-th optimal 
with respect to the class A™ nt U A^ pt as claimed in Theorem 1. 

Proof of Proposition 3 in Section 4.1 

Proof. If the problem is polynomially tractable then there exist constants C, p > and q > 
such that for all d G N and e G (0, 1] 

n(e, d) = n(e, d; F d ) < C ■ e- p ■ d q . 
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Here, e^ lt = e(0, d) = ^/\d,i > denotes the initial error of T^. Since e(n, d) = ^/\d, n +i it 
is n(e,d) = #{? G N| \d,i > £ 2 } and therefore Xd,n(e,d)+i < £ 2 - The non-increasing ordering 
of (\ d ,i)ien implies 

^d,[Cs-Pdi\+l < £ 2 - 

If we set i = [C ■ e~~ p ■ d q \ +1 and vary e £ (0, 1] then i takes the values \C ■ d q \ + 1, 
[C ■ d q \ + 2, and so forth. On the other hand, we have i < Ce~ p d q + 1, which is equivalent 
to e 2 < (Cd q /(i - I)f/P if % > 2. Thus, 

/ Cd q \ 2/p 

\ d ,i < X d ,n(e,d)+1 <Z 2 < ( — J ) fOT a11 * > maX { 2 > LC • d q \ + 1} . 

Choosing r > and f(d) = [(1 + C) ■ d q ] > max {2, [C ■ d q \ + 1} we conclude 

i=/(d) i=/(d) V 7 »=/(d)-l 

where £ denotes the Riemann zeta function. In other words, if r > p/2 > then 

l/r 



\i=/(d) / V y y 



l/r 

< oo for every del 



Setting r = 2g/p proves the assertion, as well as the claimed bound on C T . 
Conversely, assume now that (18) holds with 



/(d) = C-(rmn{ef\l})~ P -d q 
That is, for some r > and r > we have 



where C > and p, g > 0. 



l/r 



< C 2 = SUp - | > j \ T d i | < OC, 

d&i 



For n > /(d), the ordering of (A d>i ) ieN implies XT=/(d) Ki,i > Ki,n ' ( n - f( d ) + !)■ Hence, 
\d, n ■ (n - f(d) + 1) 1/T < ( £ ^ ( E ^ ^ <T, 



v<=/(<0 / V=/(«9 
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or, respectively, \ d , n+1 <C 2 d r ■ ((n + 1) - f(d) + l)" 1/r , for all n > f(d) - 1. Note that for 
e G (0, min {4 nit , 1}] we have C 2 d r ■ ((n + 1) - f(d) + l)" 1 ^ < £ 2 if and only if 



n > n* 



c 2 <r 

e 2 



+ f(d) - 2. 



In particular, it is A dj „ + i < e 2 at least for n > max {n*, f(d) — 1}. Therefore, for every dGN 
and for all e G (0, min {4 nit , l}] it is 



n(e, d- F d ) < max {n*, /(d) - 1} < /(d) - 1 + 



C 2 d rX 



5 2 



< C • (min {e™*, 1 }) P ■ d q + C 2 r 5" 2r d rr 

< ((7 _|_ (7J) . £ -max{p,2r} _ ^max{g,rr} 

Hence, the problem is polynomially tractable. ■ 

An explicit proof of Lemma 4 in Section 4.2 

Proof. Step 1. By induction on s we first show for every fixed m G N 

s 

J2 ^ k = ^ + 1^™' J2 ^ W f ° ra11 SGN - ( 41 ) 

km 8 , i=i jWen 1 , 

m+l<Jl '<...<Ji 

Easy calculus shows that this holds at least for the initial step s — 1. Therefore, assume (41) 
to be true for some s G N. Then 

oo oo 

Ms+l,fc = ^ Mfci ^ Ms,fc = Mm ^2 ^ s ' h + ^fci ^2 ^ s ' h 

k£W +1 , k x =m heN s , heW , fci=m+l fteN s , 

m<fci<...<fc 3 +i ki<hi<...<h s m<hi<...<h a k\<h\<...<h s 

= / f m ^ Ms,/i + ^ Ms+l,fc 

m<h 1 <...<h s m+l<k 1 <...<k s+1 

Now, by inserting the induction hypothesis for the first sum and renaming k to 
remaining sum, we conclude 

s 

^2 Hs+l,k = Mm 1 + ^2 Mm ^2 ^ij(') + 5^ Ms+lj(s+D- 

m<fci<...<k.+i m+l<jf ) <...<jf m+l<jf +1) <...<il; + 1 1) 



48 



Hence, (41) also holds for s + 1 and the induction is complete. 

Step 2. Here we prove (24) via another induction on V £ Nq. Therefore, let d e N be 
fixed arbitrarily. The initial step, V — 0, corresponds to (41) for s = d and m — 1. Thus, 
assume (24) to be true for some fixed V G Nq. Then it is 



/ 



i<fci<...<fc d 



l=i j(- L )eN- L , 

V+2<j[ L) <...<? j L) 
/ 



l + y + J^-* 



L=l 



V 



AV+2 



V 



1=1 jWen 1 , 

(V+2)+l<j[ l) <...<jj l) 



J 



J 



using (41) for s — L and m = V + 2. Now we estimate 1 + V by d(l + V), take advantage 
of the no n- increasing ordering of (// m ) me N and extend the inner sum from L to d in order to 
obtain 



/ 



^d,k < (4 d 

k&i d , 
l<kx<...<k d 



d iV+1 



\ 



i+(v+i)+j2»i l E 



Hu(i) 



\ 



i=i j(0 e N ( , 

(V+l)+2<j<°<...<j< 



Since this estimate corresponds to (24) for V + 1 the claim is proven. 



Proof of Proposition 8 in Section 4.4 

Proof. Note that due to lim™-^ \ m = the quantity id is well defined, because (ai)i G N is 
a non- increasing sequence which tends to zero for i tending to infinity, and id(S 2 ) > 1 for 
5 < e™ 1 *. Furthermore, we have 

i^S 2 ) = #{m G N | \ m >5 2 } + 1 = n cnt (5, 1) + 1 = n asy (5, 1) + 1 

and if d > 2 we can rewrite ^ to obtain 

i d (S 2 ) = min j i e N | A i+ i ■ . . . • \ i+ d-i < ^ 
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Hence, for every k = (ki, . . . , fc<2-i) G V^-i with k\ > id(8 2 ) it is 

Ad_i,fc = A fcl • . . . • A fcd _ x < A id(5 2 )+1 • . . . • A, d (52)+d-i < t S 2 

or, equivalent ly, 

k G V d _i | i < fci and A d _i jfc > -^5 2 1=0 for all i > i rf (5 2 ). 

Aj J 

This leads to the disjoint decomposition of 

{j G V d | A dJ > 5 2 } = | j = G N x Vd_i|i < h and A d _ life > ^-5 2 | 

= M < (i, fc) | fc G Vrf_i such that i <k\ and A rf _i fc > — S 2 > . 
i=i ^ Al J 



Therefore, the information complexity of the d-variate problem is given by 

n^(e,d) = #{j'G V d \X d , J >e 2 }= £ # |fc G V d _! \i < h and A d _ life > -e 2 j 

= 52 # \ k e V ^-i I z i < and A <i-i,fe > t e 2 \ . 



h=2 

_2 



Obviously, for fixed l\ G {2, . . . , )}, we can repeat this procedure and obtain 
#{j G Vd_i | Zi < ji and Ad_ij > 5 2 } 

i d -i(5 2 ) 



V # (fc G V d _ 2 \l 2 <h and A d _ 2jfc > y^5 2 ) 



/ 2 =/i+i 

if d > 2 and 5 2 = e 2 /\ h _ 1 . Note that £ < ejf'* implies ^_i(5 2 ) > h + 1 such that {/i + 
1, . . . , id-i(5) 2 } 7^ 0. Iterating the argument we get 

n asy (e,d) 

i d (e 2 )id-i(e 2 /\ h -i) i2(e 2 /[X h - 1 :.,X ld _ 2 - 1 ]) 

= J2 Yl ••• Yl # <| ^ e Vx | Z d _! < fci and Ai ifc > - -<5 2 

h=2 l 2 =h + l i d _i=J d _ 2 +l 
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with 5 2 = e 2 /[Ajj_i ■ . . . • Aj d _ 2 _i]. It remains to calculate the cardinality of the last set. Of 
course, we have 

k e Vi | Z d _i < fci and Ai )fc > — — — ^ 2 

fc G N | l d -i < k and A fc > — — — 5 2 

A/d-i-i 

fc e N | A fe > v^5 2 ) \ ( k e {1, . . . , Z d _ x - 1} | A fc > v^5 2 



The first of these sets in the last line contains exactly n ent (<5/ A yA/ d _ 1 _i, 1) elements. On the 
other hand, if k < la-i < ii{& 2 ) then 

AfcA; d _ 1 _i > \ 2 {8?)\ 2 (jp)-i > 5 2 , 

where the last inequality holds due to the definition of i2{8 2 )- Therefore, the last set coincides 
with {1, . . . ,ld-i — 1} and its cardinality is equal to l d -\ — 1- Furthermore, note that the 
estimate also shows that n ent (S/ ■^/Aj d _ 1 _i, 1) is at least equal to ld-i- Thus, 

# \k e Vi | l d -i < fcx and A 1)Jfc > = ^ (Vy^-i-i, l ) - k-i + 1 > 1 

and the proof is complete. ■ 

Proof of Proposition 9 in Section 4.5 

One possibility to prove the second point of Proposition 9 is to apply proposition 3 to a 
scaled problem {T d } such that Wd = T\T d possesses the eigenvalues A^j = A^j/A^i for 
i G N. Then the initial error of Td equals 1 such that / in Proposition 3 does not depend 
on p. That is, we can choose f(d) = \C d q ] + 1 for some q > in both the assertions. In 
order to see why we even can take f(d) = 1 in the first point and for the sake of completeness 
we also add a direct proof for this proposition. 

Proof. If {Td} is polynomially tractable with respect to the normalized error criterion then 
there exist constants C,p > and q > such that for all d e N and e' G (0, 1] 

n(e' ■ ef\ d) = n(e' ■ ef\ d; F d ) < C ■ (e'y p ■ d q . 



As before the quantity e™ 1 * = y/Xd,i > denotes the initial error of T d and e' is the (mul- 
tiplicative) improvement of it. Since e(n,d) = ^/\ d , n +i it is n(e,d) = € N | \ d ,i > e 2 } 
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where e = e' -e 1 ^. Therefore, A djn ( e /. e init jd ) +1 < (e' f ■ \ djl . Hence, the non-increasing ordering 
of (Ad,i)igN implies in this setting 

X d,[C(e')- p d"\+l ~ ( £ ') 2 ' Xd ^- 

If we set i = [C (e')~ p d q \ + 1 and vary e' G (0, 1] then i takes the values [Cd q \ + 1, [Cd q \ + 2 
and so on. Again we have 1 < % < C (e') p d q + 1 on the other hand, which is equivalent to 
{e'f < (Cd q /(i - l)f/P if % > 2. Thus, 



A(i,i < A d)fl ( e /. e init )d ) + 1 



/ nM \ 2 /p 

< (ef • A d) i < f — - J • A dil for all i > max {2, \Cd q \ + 1} 



Choosing r > and f*(d) = \(1 + C) d q ] > max {2, [Cd q \ + 1} we conclude here 
where C, again is the Riemann zeta function. On the other hand, it is obvious that 

r(d)-i 



E J < r (d) - 1 < (i + c) ^ 



^ \Ad,i 

because A di j < A dj i for all ieN. Therefore, if r > p/2, 



oo 



for all d G N. This proves the assertion setting r > 2q/p. 

The proof of the second point again works like for Proposition 3. Assume that (37) holds 
with f(d) = \C d q ~\ , where C > and q > 0. That is, for some r > and r > we have 



d ^ d V^) V di1 



1/T 

< OO. 



Since (A d) j)j S N is assumed to be non- increasing the same also holds for the rescaled sequence 
(A d ,i/A d) i) i6N such that T / "=f(d)( X d,i/^d,i) T > (K,n/^d,i) T • (n- f(d) + 1) for ra > /(d). Hence, 

... 1 / T 

A<2,n 



Ad,i 



n-/(rf) + i)^< ^f; <c T <r, 
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or, respectively, \d, n +i < C T d r ■ ((n + 1) — f(d) + 1) 1 I T ■ for all n > /(d) — 1. As before 
we have C T d r ■ ((n+ 1) - /(d) + l)" 1 ^ < ( £ ') 2 , for e' E (0,1], if and only if 



n = 



C T d r 
W 



+ /(d) - 2. 



In particular, A^ n+ i < (e') 2 • A d l at least for n > max {n*, /(d) — 1}. Therefore, we conclude 
in this setting for all e' G (0, 1] and every d e N 



n( £ ' • er , d; F d ) < max {n*, /(d) - 1} < /(d) - 1 
< (C + C T T ) ■ (e'Y 2T ■ d max ^' rT \ 



^) <Cd 9 + C T r (^)" 2r d r 



Hence, the problem is polynomially tractable. Furthermore, strong polynomial tractability 
holds if r = q — 0. ■ 
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